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ABSTRACT 

The problem of designing robust con 'ol systems for the safe landing 
of VTOL aircrafts on small ships is addressed for the lateral motions. 

A precise ship model is derived, using hydrodynamic data for the DD963 
destroyer. The issues of estimation and prediction of ship motions are 
discussed, using optimal linear estimation techniques. The roll motion 
is the most important of the lateral motions, and it is found that it 
can be predicted for up to 10 seconds in perfect conditions. 

The automatic landing of the VTOL aircraft is considered, and a 
lateral controller, defined us a ship motion tracker, is designed, using 
optimal control techniques. The tradeoffs betweei, the tracking errors and 
the control autliority are obtained. The important couplings between the 
lateral motions and controls are demonstrated, and it is shown that the 
adverse couplings between the sway and the roll motion at the landing jiad 
are significant constraints in the tracking of the lateral ship motions. 

The robustness of ti;^ control system, including the optimal estimator, 
is studied, using the singular values analysis. Thiraugh a robustification 
procedure, a robust control system is obtained, and the usefulness of the 
singular values to define stability margins that take into account general 
types of unstructured modelling errors is demonstrated. The minimal desta- 
bilizing perturbations indicated by the singclar values analysis are 
interpreted and related to the multivariable Nyquist diagrams. 
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CHAPTER X 


INTRODUCTION 


X.l Introduction and Motivation 

The landing of small VTOL aircraft on destroyers is an extremely 
challenging problem if it is to be realized under high sea state conditions 
and zero visibility. Without special aids, this task is almost impossible 
for a human pilot. 

There arc basically two possible strategies in the solution of this 
problem. The first is to leave to the pilot the complete control of the 
aircraft, but help him with advanced displays. These give him information 
about the aircraft position and attitude, as well as those of the ship 
(and possibly some prediction of the ship motions) . They may also indicate 
some desirable flight path (flight director) . Advanced controls may be 
provided to partly relieve the pilot from the high load of controlling the 
VTOL aircraft. 

The second strategy is to leave the task of landing the aircraft 
completely to an automatic controller. The role of the pilot is then to 
supervise the correct landing of the aircraft. This wou].d allow him to 
take care of other tasks he might not have been able to carry out other- 
wise. 

Note that both strategies could be mixed. For example, the tracking 
of the lateral ship motions may be left to an automatic controller/tracker, 
while the task of vertical tracking and landing would be left to the pilot. 
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possibly with the help of some display indicatinq him the present and 
future ship vertical position. 

In this thesis, the emphasis is focused on the design of an automatic 
controller. A previous study (Ij has addressed the problem of the longitu- 
dinal motions, i.e. the motions in the vertical plane. The most significant 
ship motion in this plane is the vertical motion, called heave. The pitch 
motion is not very significant, except for the heave motion it induces at 
the landing pad (which is significantly behind the ship center of rotation). 
The present study addresses the ship ijotion tracking problem for the lateral 
case. Then, the most significant motion is the ship roll motion, which can 
be very large. The lateral translation motion, called sway, is also impor- 
tant, especially due to the large sway component induced by the roll at the 
landing pad (located above the ship center of rotation) . 

The challenge of the tracking of the ship motions by the VTOL lies in 
the strong limitations of the control authority available, in the high 
level of the perturbations (wind disturl>ances, ground effects, ship airwake), 
in the strong couplings present in the system, and in the need for a highly 
robust control system. 

Usually, studies of this problem use loop-by-loop control system 
designs, using classical control theories [ 2] , [3] , [4] , [5] , [6] . In this 
case, the controller ignores the internal couplings of the system. Similarly, 
the issue of robustness is often addressed on a loop-by-loop basis, but 
almost never in a real multivariable sense (although individual loop sta- 
bility margins may not represent at all the overall system stability 
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mai'gina) . The design process used in this thesis does not sufie} fiom 
such iimitationa. 

The iimitations on the available control autliority justify the use 
of some optimization criterion, and of related modern control theories 
(IjQ/LQG) . These methods have the advantage of naturally handling mult i- 
variable systems, and of recognizing the coupling present in such systems. 
Some recent results in the analysis of the robustness of multivariable 
systems (and its improvement for LQG based designs) are also important 
tools in the c osign of control systems operating under critical condi- 
tions . 

The purpose of this thesis is not to produce an engineering design. 
Nor does it provide new theoretical results. It illustratet: liow modern 
control theories and related recent results can be used to design a 
control system for such an advanced application, and evaluate the 
controller performance and robustness. This work also analyzes the physi- 
cal constraints of the tracking process of the lateral ship motions. 

These constraints are independent of the control system design methodo- 
logy adopted. The requirements and physical limitations related to the 
VTOL landing pr-oblem are studied. 

Although this thesis mainly details the design of an automatic 
controller, the accuracy achievable in the prediction of ship motions is 
also assessed, as it is a key element in any piloted VTOL landing. 


1.2 ContributionB of the Research 


The main contributions of the research are: 

- the derivation of an accurate ship model that retains the stochastic 
nature of the ship motions, and the couplings amongst them 

- the analysis of an optimal predictor of the ship motions for 
apiJiications in piloted landings, and the assessment of lower bounds on 
the prediction errors 

- the design of an optimal controller/tracker for applications in 
automatic landings, the indication of the tradeoffs between tracking errors 
and control authority, and the analysis of the important couplings and 
physical constraints related to the tracking of the lateral ship motions 

- the demonstration of the use of the singular values analysis, and 
the robustification procedure, to obtain a robust control system. 

1 . J Thesis Overview 

Chapter 2 indicates the general problems of the landing of VTOL 
aircraft, and introduces the methodology used to design the control system, 
as well as some aspects specific to this particular application. 

Chapter 3 details the derivation of the shiis model from hydrodynamic 
data, and indicates the important characteristics of the model. 

Chapter 4 is an analysis of the issues of estimation and prediction of 
the ship motion, using the model derived in chapter 3. 

Chapter 5 summarizes the aircraft model obtained from [ 1 ] , and gives 
a brief description of this model, and of the important couplings present 
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between the motions and controls. 


Chapter 6 details the design of an optimal control system based on 
LQ theory to track the ship motions. Important characteristics such as 
root-locus, step responses, and tracking errors versus control authority 
are analyzed, with their relation to the couplings amongst the motions 
and controls. 

Chapter 7 shows the design of an optimal estimator for the aircraft 
motions, and indicates the degradation of performance due to the presence 
of noise in the sensor measurements. 

Chapter 8 addresses the important issues of robustness of the control 
system to modelling errors, and demonstrate the usefulness of the singular 
values analysis and of the robustif ication procedure. 

Chapter 9 concludes with some general comments and suggestions for 
further research. 
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CHAPTER 2 


OVERVIEW : OBJECTIVES AND METHODOLOGY 


2 . I Introduction 

The landing of VTOL aircrafts on small platforms is a delicate 
operation that interests civilians (oil platforms) as well as military 
(destroyers) . In this research, we consider a moving platform, specifi- 
cally the landing pad of a destroyer, type DD-963, in sea state 5. Such 
sea state corresponds to waves of heights around 10ft, and winds around 
20 kts. Until now, this is still a goal, and it justifies studies as 
this to assess the navigation systems, the aircraft performance, and 
the control system required to perform such an operation. 

A previous study by McMuldroch [1] has addressed the VTOL landing 
problem for longitudinal motions. These are called heave (vertical 
motion), surge (fore and aft motion), and pitch. The aircraft considered 
in this work (as in the present work) is the Lift/Pan Cruise Research 
Technology Aircraft (RTA) . More details about this aircraft can be found 
in I l) and in chapter 4. The specific aircraft studied is not really of 
great importance, as most VTOL's have the Scune kind of limitations and 
possibilities. The main characteristic of an RTA-type of VTOL aircraft 
- which makes it different from a helicopter for example - is the possi- 
bility of deflecting the engine thrusts to produce translation motions 
without rotating the aircraft. In other words, this kind of VTOL has as 
many controls as it has degrees of freedom, so that, to the limit, 


- 19 - 



perfect tracking of the landing pad motion can be acheived. 

In the present work, we use the model derived in 111 for the RTA 
aircraft. Our goal is to complete McMuldroch's work with a study of the 
controls needed for the lateral motions. The lateral motions are called 
sway (lateral translation), roll and yaw. The longitudinal and lateral 
tiKitions can be decoupled to a good approximation for the ship motions. 

This does not mean that they are independent; in fact, they are strongly 
correlated because they are generated by the same wave. However, except 
for this commc n excitation force the dynamic equations for longitudinal 
and lateral motions can be decoupled to first order and can be analyzed 
separately. The decoupling of longitudinal and lateral motions for the 
aircraft is a little less obvious, due mainly to gyroscopic cross-cou- 
pling terms [7] . These terms can be quite important, due to the large 
size of the engines of a VTOL aircraft. To a first approximation, however, 
they can be neglected, so that longitudinal and lateral motions can be 
studied separately. 

In fact, the issues for the longitudinal and for the lateral case 
are quite different. For the longitudinal case, the important motion is 
heave. Pitch is quite small and surge is negligible. Moreover, they can 
be controlled easily. Heave is more critical, as it requires an increase 
in the overall engine power (or thrust to weight ratio). This is a control 
which is strongly limited in amplitude, and also in its speed of response. 
The limitations are so strong that one may look for an end-point 
controller, instead of a tracking controller. 
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In the lateral case, two of the three motions are about as important: 
they are the away and the roll motion. The yaw motion is almost insignifi- 
cant, and can be controlled very easily. In many situations, the aircraft 
roll may bo left uncontrolled, or, more precisely, stabilized to zero 
(so that it does not track the shlp% roll) . In this case, the landing gear 
will simply damp the shock at the landing. However, ship roll can be as 
large as 30 degrees peak to peak (even more in the case of decaying seas), 
so that this strategy is not always satisfactory. Then, both ship roll and 
svray motions have to be tracked. 

The limitations in control authority are less severe than in the 
longitudinal case. In particular, a lateral side force can be produced by a 
deflection of the thrusts, and a roll moment by an exchange of power 
between the engines. These controls can be generated very quickly, so 
that there is no significant problem with bandwidth in the lateral case. 

The problem is probably more in the adverse couplings present in the 
lai'.«?.ral motions and in the robustness of the closed-loop system. These 
tvi \;,pects will be examined in more detail later. The consequence of 
thes ■ remarks is that the problems are slightly different for the longitu- 
dinal and lateral cases. The lateral controller can be a tracking 
controller, while the landing itself occurs at a moment decided by the 
longitudinal controller. The responsibility of the lateral controller is 
then to minimize the tracking errors at impact. 

As the research on the longitudinal motions had shown that precise 
ship modelling was essential to obtain meaningful results, a large part 
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uf this thesis is devoted to accurate ship motion modelling, estimation 
and prediction. The prediction part of the study is irrelevant to the 
automatic controller design part, but has applications in piloted landings, 
when some indication has to be given to the pilot about where the ship 
deck is now and will be in the near future. 

2 . 2 Control System Design Methodology 

2.2.1 Inj^rod^c^io^n_ 

Differen< methodologies have been proposed recently for the design 
of multivariable control systems. Whether applied in the time domain or 
in the frequency domain, each metfod has its own advantages and limita- 
tions. 

For aerospace applications, the linear-quadratic-gaussian (LQ-LQG) 
methodology has shown successful, especially because; 

- the limited control authority available makes the optimization 
imperative 

- the number of states is small (with a good approximation) 

- the equations of motion (and of the system in general) are quite 
well known, and the state-space description is natural 

- the systems are often unstable, and strongly coupled. 

2.2.2 bP_Methodol^ 03 Y_ 

We summarize in this section the LQ methodology to specify the 
notation used subsequently. More detailed descriptions can be found in (8] 
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Given the system: 

^ Ax Bu + ^ (2.1) 

x(0) ■= ^ (2.2) 

whore : 

X is the state vector (nxl) 

X is the Initial state vector (nxl) 

~o 

u is the control vector (mxl) 

C is a white Gaussian noise vector (nxl) (driving noise) 
of spectral intensity matrix (nxn) E 

A is a matrix (nxn) 

B is a matrix (mxm) 

One wants to minimize the quadratic optimization criterion: 

J « lim E { ^ (x^ Q -f R u) dt } (2.3) 

T-H» 0 

where: 

Q is the positive semi-definite matrix (nxn) of the state weights 
R is thf? :x>sitive definite matrix (mxm) of the control weights 
The pair (A,B) is assumed to be stabilizeable and the pair (A^Q*^) 
detectible (two conditions easily satisfied). 

The solution is a time-varying control law u(t) which, for large t 
(far from the origin of time), becomes a linear, time-invariant. 
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state feedback: 


y. “ " G X (2.4) 

where G is an (mxn) matrix given by: 

G •• r“^d’’k (2.5) 

and K is the (nxn) positive matrix solution of the matrix algebraic 
Riccatti equation: 

KA -t a'^^K - KBr'^b'^K + Q » 0 (2.6) 

2.2.3 I^G Methodology 

The LQG methodology is the extension of the LQ methodology when the 
state X is not available for measurement but, instead, we have a measure- 
ment vector y, which is a noisy linear combination of the state 

y “ Cx 4- (2.7) 

where; 

y is a measurement vector (pxl) 

C is a matrix (pxn) 

£ is a white Gaussian noise vector (pxl) (measurement noise) 
of spectral intensity matrix (pxp) 0 

The optimization problem is the same and, from the separation principle, 

it is known that the solution is the cascade of an optimal estimator 

(Kalman filter) providing an estimate k of the state vector x, and of 
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the same control law as before, acting on ^ instead of xr 

H. * " G B. (2.0) 

The Kalman filter structure is: 

^ 4- Bu + (2.9) 

where H is an (nxp) matrix given by: 

H » P c'^ 0’^ (2.10) 

and P is an (nxn) matrix, positive solution of the matrix algebraic 
Riccattj equation: 

Pa'^ + AP 4- 5 - PC'^0"^CP = 0 (2.11) 

The solution of Riccatti equations is now done routinely by the use of 
a modern control system design package [10] . 

2,2.4 hp_C 2 .*^trol^e£ and S^h^_T£ack£ng_ 

In the VTOL landing problem, we are interested in the tracking, 
by the aircraft, of certain ship states (ship motions and velocities). 

We have a situation in which we wj to track an uncontrollable system 
(the ship deck) by a controllable system (the VTOL). 

For the ship, we have the following state-space description: 

Bs=AgXg4£g (2.12) 
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(2.1J) 




vAiere are the ship states that we want to track. 

“k> 


For the aircraft: 


X, “ + B,u 4- £, 

A-A A— 


YL. ~ 

A-^ 


where w. are the aircraft states corresponding to 

—A ~b 

The general system equations are: 


X = Ax + Bu + £ 


(2.14) 

(2,]5) 


(2.16) 


where: 


X = 



A 



(2.17) 


(2.18) 


B 


0 



i 



(2.19) 


( 2 . 20 ) 
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We want to minimize the tracking errors and the control authority, i.e 


J - + u^Ru) dt) 


o 


,00 T 1' 

E(J ()£ Qx + ^ Ria) dt) 


where: 


Q = 


T 

WO W 

s'^w s 

-wrc 

Aw S 


T 

-W„Q W„ 
S^w A 

T 

WnQ W, 

Aw A 


The optimal LQG controller is: 


( 2 . 21 ) 

( 2 . 22 ) 


(2.23) 


- G X = - 


‘°s ' 


A 


- G^x^ - G, X- 
S-S A -A 


(2.24) 


Since the ship is uncontrollable, the ship states and the gain 
matrix G appear as a feedforward path, while the feedback path is 
coming from the aircraft states only, through the gain matrix G . 

ri 

As indicated in the next paragraph, the estimation problems are also 
decoupled, so that the controller structure is that of Fig. 2.1. 


L 



Figure 2.1 : VTOL Controller Structure 
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Ail important remark is that the gain matrix G is the same than the 

A 

matrix obtained by solving the optimization problem.* 

+ u^Ru) dt) (2.25) 

O 

The gain matrix G is then independent of the ship model. 

'To see this, we assume that the solution of the Riccatti equation 
corresponding to (2.21) is ; 


K 


'■'ll ■'12 

■'L *=22 

L 


Then s 


-1 T 
R B K 


(2.26) 


(2.27) 


-1 T 

- I" Vl2 


■'■*❖ 22 ’ 


(2.28) 


and: 


=A " ■'■'‘’^22 < = 

depends only on ^^ 22 ’ 

The Riccatti equation is (2.6), and can be partitioned in terms of the 
ship and the aircraft parts: 


T -1 T T T 

K, + A^K, , - K, .,B,R B,K, ^ + WO = 0 
11 S S 11 12 A A 12 S^W S 


^12^ ^S*^12 “ ^12®A^ ®A^22 ~ '^sVa ° 


(2.30) 

(2.31) 
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(2.32) 


*^ 22 \ ^ \*^22 


- '^ 22 V ‘'“>22 “ “ 


The equation giving K is independent of the others and is actually 
Die same as the one corresponding to the quadratic criterion (2.25). The 
same is true for the gain matrix G^. This means that the feedback gain 
matrix G^ is independent of the dynamics of the signal that is tracked, 
namely the ship in this case. The ship only influences the feedforward gain 
matrix. In other words, the closed-loop dynamic behavior of the system is 
only dependent on the Q and R matrices, and is the same whether we want to 
track the ship or simply to stabilize the aircraft (track a zero reference 
signal). For this reason, the closed-loop behavior fo the system can be 
studied independently from the ship model (in this, we include the optimal 
root-locus, the step responses, and the robustness characteristics of the 
system) . Note liowever that the ship model influences the system performance 
(tracking errors) which, in turn, may influence the specific choice of the 
Q and R matrices. 

2.3 State Estimation and Implementation Issues 

The accuracy of the measurements of the ship and aircraft motions will 
be a determining factor of the performance attainable in the landing of VTOL 
aircraft in high sea states and poor visibility conditions. It is expected 
that, in a practical realization, the measurement process will involve: 

- accelerometers and gyroscopes aboard the aircraft 

- a combined microwave landing system / distance measuring equipment 
(MLS/DME) giving measurements of the relative position of the ship and the 
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the aircraft. 


TheoreticaJ ly, the estimation of the aircraft and of the ship motions 
are coupled, in the sense that the relative position measurement would 
influence the estimation of the ship motions if a Kalman filter is used 
for the overall system. In practice, the estimation of the ship motion, 
and the estimation of the aircraft motion can be reasonably decoupled. 

A justification for this is that the instruments aboard the ship will be 
of somewhat greater precision than those aboard the aircraft, and the 
estimation errors obtained will be smaller than the errors on the aircraft 
motion. The estimation part of the LQG algorithm can then be separated in 
an estimation of the ship states by a Kalman filter, and a similar but 
independent estimation process for the aircraft. 

Practically, the estimation of the ship states is done aboard tiie 
ship, and the results are data-linked to the aircraft. With the combined 
use of the accelerometers, gyroscopes, and MLS/DME measurements, a 
computer aboard the aircraft estimates the aircraft-ship relative positions 
and velocities, and through the use of a Kalman filter, obtains optimal 
estimates of the aircraft states. The control law is then easily obtained 
aboard the aircraft. 

As the computing capabilities aboard the ship can be more powerful, a 
high order model can be used. In our case, it is a 16 states model (see 
chapter 3). On the aircraft, the computing capabilities are more limited, 
so that a low order model is desirable. The aircraft model used here is a 
6 states model (see chapter 5). The importance of robustness becomes 
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critical, however. While the approximations on the ship model have conse- 
quences only on the performance of the controller, those on the aircraft 
model come in the feedback structure, and influence the robustness as well 
as the performance of the controller. 

2.4 Summary 

In this chapter, we discussed the general aspects of the VTQL landing 
problem. The decoupling between the longitudinal and lateral motions was 
explained, and it was indicated that the controller for the lateral motions 
can be conceived as a ship motion tracker. 

The LQ/LQG control system design methodology was introduced, and its 
use for this application was justified in view of the limited control 
authority available. It was demonstrated that the application of this 
methodology to the VTOL landing problem leads to a decoupling of the ship 
and of the aircraft effects, so that the ship motion and the ship dynamics 
only influence the feedforward structure of the control system, while those 
of the VTOL impact the feedback structure. 

Finally, the issues of estimating the ship and aircraft states were 
briefly addressed, and it was shown that the problems of estimation of 
ship and aircraft motions could be reasonably decoupled. 
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CHAPTER 3 


SHIP MODEL DERIVATION 


3.1 Introduction 

The importance of accurate ship modelling in VTOL landing has been 
indicated by McMuldroch in [1] . Ship motions have relatively narrow band 
power spectra (between 0.2 and 2 rad/s) that require high order models 
to be represented accurately. A good ship model, that is often used in 
studies of VTCL landings, is available in [11] . The motions are modelled 
as sums of 6 to 32 sinusoids, with random phases. For simulations (espe- 
cially for flight simulators), such a representation is adequate. 

However, it fails to represent the random nature of ship motions for 
longer periods of time, and it is not appropriate for estimation of ship 
motions from noisy measurements, and for ship motion prediction. 

In this chapter, we derive a ship model that retains the stochastic: 
nature of the ship motions and the important couplings among the various 
motions. The model is derived in state-space form, so that the powerfu] 
techniques of linear estimation in the time domain can be applied. The 
equations are obtained from hydrodynamic considerations, that lead to 
linear differential ecjuations with frequency dependent coefficients, and 
infinite dimensional transfer functions. Finite dimensional approximations 
are considered, and the model finally includes 16 states for the lateral 
motions. The following sections indicate the structure of the model and 
the approximations made. More details can be found in [12] , [13] and [14] . 
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3.2 Ship Mcxiel Structure 


The lateral ship motions, and the sign conventions used here, are 
shown on Fig. 3.1. The origin of the reference system is located in the 
plane of symmetry of the ship, at the level of the waterline, and at the 
middle of the ship. This point is close to the center of gravity and has 
the advantage of being independent of the ship load conditions. 

The ship motions are assumed to be small, in order to derive linear 
equations of motion. This assumption is justified by the fact that waves 
have limited wave to length ratio (at most 1/7, and usually much smaller), 
since waves having higher heights break and loose their energy. As a 
result, the major part of the force is linear and can be obtained by a 
first order perturbation expansion of the non-linear fluid equation. 

The wave spectrum is typically contained in the 0i2-2 rad/s range. Given 
the large mass of the vessel, the resulting motions, within this Irequency 
range, are of the order of a few feet and a few degrees, so that the 
linearity assumption can be justified. Roll motion requires more attention 
however. Due to the slender form of the ship, the roll motion may become 
large and then, the non-linear damping is predominant. 

The ship model is basically divided into tliree parts. The first part 
represents the incoming waves (sea model) which are described by the wave 
elevation at a reference point located amidships. The wave elevation is 
known to be a stochastic process, defined by a relatively narrow band 
power spectrum to which various approximations have been proposed. 
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The second and third parts of the model put together, represent the 
runponse of the ship to tJie waves. Using the: hyiKJthesis of linearity, this 
response can be represented by a transfer matrix from the wave height to 
the ship motions. The second part of the model represents the dynamics of 
the forces generated by the incoming waves, while the third part of the 
model represent the rigid body dynamics as well as the dynamics induced by 
the ship motions on the water motion resulting in additional fluid forces. 
These forces are : 

- the inertia force caused by the acceleration of the fluid particles 
displaced by the ship (added mass) 

- the damping force caused by the loss of energy carried away by waves 
generated by the ship motions 

- the hydrostatic force (spring constant). 

Por a purely sinusoxdal wave, the equations of motion lead to the follo- 
wing differential equation s 


(M+A. )x. + bJc. + C. X. = F 
h — h h-h h-h — 


( 3 . 1 ) 


where : 

- M is the mass matrix of the ship, including mass terms, products of 
inertia and coupling terms due to the difference between the center of the 
axes and the center of gravity 

- A. is the added mass 

h 

- is the vector of the ship motions (sway, roll and yaw) 

- is the damping term 
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- is the hydrostatic term 

- F is the vector o£ the forces and moments generated by the incoming 
wave 

The terms F, A^, and depend on the frequency of the incomj.ng wave, so 
that the equation of motion is strictly valid for monochromatic waves. 

For an irregular sea elevation, equation (3.1) becomes an integro-dif fe- 
rential equation, hh additional difficulty is due to tte fact that wave 
forces and moments are obtained by integrating over the ship hull the 
space-varying pressure forces, so that their magnitude and phase constitut 
the transfer function of an infinite dimensional system. Obtaining a finite 
dimensional model of reasonable size and complexity definitely requires 
several simplifying approximations that will be indicated in the subsequent 
paragraphs. 

The geometric and mass properties of the DD-963 were analyzed by the 
M.I.T. Ocean Engineering Department Seakeeping Program [15] . The hydrody- 

namic coefficients (M,A. ,B, ard C.) and forces were first obtained, and 

h h h 

subsequently the overall ship model was derived. 

The parameters of the model are : the speed of the vessel, the direc- 
tion of the waves, the significant wave height, and the modal frequency of 
the wave spectrum. 

3.3 Sea Modelling 

The sea waves are generated by the wind, except for very rare cases 
(seismi waves). The high frequencies of the wind gusts create wavelets on 
the surface of the sea, while the steady-state condition of the sea develops 
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Blowly t-'\rough a nonilnoar interaction mechanlBm creating waves whose 
phase velocity is close to the wind speed. Since the process starts with 
high frequencies, a young storm will contain a peak at high frequency. As 
soon as the wind stops blowing, the water viscosity dissipates the high 
frequency waves so that the so-called "swell" forms, which consists of 
long waves (low frequency) which travel away from the storm that origi- 
nates them. For this reason, swell can be found together with another 
local sitorm, in which case the wave height spectrum contains two peaks. 

The intensity of the storm can be described in various ways, of 
which the best is probably the significant wave height H, defined as the 
statistical average of the 1/3 highest waveheights. At any point, the wave 
elevation is a stochastic process described by its ixswer spectrum. This 
power spectrum (if it is single-peaked) is a function primarily of two 


parameters, H and to^. The Bretschneider spectrum is defined as : 


S(ai) = 


0 )^ 


U) 


exp (-1.25 (0)^/w)^) 


(3.2) 


It was found to fit reasonably well in any sea location, and is strictly 
valid for unidirectional seas, with unlimited fetch, infinite depth, and 
no swell. 

As the ship moves toward the waves, the apparent frequency of the 
wave is modified, and the frequency of encounter is ; 

= 0) + k U cos(j) (3.3) 
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where : 


- U is the ship speed 

- i() is the angle between the x axis of the ship and the direct ion 
of the wave propagation 

k is the wave number. 

In deep water, the dispersion relation for the waves is s 

0)^ * kg (3.4) 


so that : 


0) - 0) 4- U COS(j) 

e g ^ 


The spectrum seen from ship coordinates is then 




S(M) 

d(i) /dw U)=f (W ) 
e •' e 


(3.5) 


(3.0) 


where ; 


- 1 + / 1 + 


4 0) 


U COS(() 


0) = f (w ) = 
e 


2 — cos(|) 


(3.7) 


To obtain a state-space representation, the wave height is represented as 
the output of a filter with rational transfer function, driven by white 


noise. 



The transfer function is selected to be : 


li (s) /ili"' 
a o 


(“) 

0) 

o 

1 2 J(§) 4 

o o 


(3.8) 


whore the values of S^, U)^, and J are dependent on the sea state and are 
qiven in Appendix A. 

The coefficient il was introduced to match the differences in defini- 
tion for the power spectrum. Here, we use : 


/<x> 

8(0)) = / B(T) e dx (3.9) 

.00 

A one-sided (positive fre<iuency) power spectrum is used in wave theory 
(and lor the Bretschneider spectrum): 


1 — l(i\T 

S( 0 )) = i f R(X) e dx 0»0 (3.10) 

~QO 

The Bretischneidor spectrum and the approximation are illustrated in 
F i.g . 3 . 2 . 

3.4 Ship Transfer Matrices 

♦ 

The ship transfer matrix from the wave height to the sway, roll and 
yaw motions can be separated in: 

- a transfer matrix from the wave height to the force and moments 
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Figure 3.2 ; Bretschneider Spectrum and Approximation 

(H=10ft, 0 ) =0.72racl/s) 
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2.00 


generated by the incoming wave 

- a transfer matrix from the force and moments to the actual motions. 
Those transfer matrices are obtained from the values of the forces^ added 
mass, damping, mass, and hydrostatic forces, using equation (3.1). For 
monochromatic waves, these were obtained by the M.I.T. Ocean Engineering 
Seakeeping program from the geometry of the DD-963 hull. Rational approx- 
imations were made in order to obtain a standard state-space representation. 

An important fact is that whatever approximations are made, they do 
not influence tha stability margins of the aircraft control system. The 
ship model only comes in the ship Kalman filter design and in the compu- 
tation of the feedforward gain matrix. The only effect of the approxima- 
tions on the control system is on the performance attained in the ship 
motion tracking, not on the system robustness. 

3.4.1 Fo£C£ 2y>}355.i£,s_ 

The transfer functions between a unit cunplitude regular wave (1ft) 
and the sway force (tons) and the roll and yaw moments (tons-feet) as 
functions of the frequency of encounter are shown in logarithmic scale 
in Pigs 3. 3, 3. 4, 3. 5 for the case U=0,(f>=90 degrees and in Figs 3. 6, 3. 7, 3. 8 
for the case U=15.5ft/s and <()=45 degrees. 

The forces were approximated by simple second-order systems, for 
exiunple : 

F. s 

P, =: i i=s,r,y (3.11) 

1 + 2 j./w. 4- (s/0).)^ for sway, roll and yaw 
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SWRY force: 



00 (RAD/s) 


Figure 3 , 3 : Sway Force and its Approximation (U'=0, 4*“^OdG9) 



Amplitude (dB) 


ROLL MOMENT 



Figure 3.4 ; Roll Moment and its Approximation (U=0, (|)=90deg) 
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Figure 3.8 : Vaw Nbment and its Approximation (Uel5.5ft/s, 
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In the case when U«15.5Ct/s, (|)=45degrces, the higher frequency behavior 
shows the infinite dimensional cliaractoristic of the force dynamics. I’or 
simplicity, this effect is not represented in the model, because it is 
beyond the area of significant wave power. 

The values of the coefficients in the transfer function, and their 
dependence upon ship speed and wave heading angle, are summarized in 
Appendix A. 


3,4.2 S^h^_Dyn^^C£ 

The other ship dynamics are related to the matrices Aj^, and 

As indicated earlier, the matrices A^^ and are dependent on the frequency 
of the incoming wave, so that the order of the differential equation: 




B, X, + C^x. = P 
Irdi h-’h ~ 


(3.12) 


is higher than 2 times the three motions. 

However, the roll motion turns out to be highly concentrated around 
the roll peak frequency. The sway and yaw power spectra are more widely 
dispersed, but, due to the concentration of the sea spectrum in frequency, 
the motions are all concentrated within a narrow frequency band. 
Consequently, the added mass and damping matrices variations with frequency 
are neglected, and the values of these matrices at the roll peak frequency 
are used here. 

Another important variation of the Aj^ and B^ matrices is due to the 
ship speed and heading angle. These produce coupling terms between the 
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raotiono that, again, are frequency dependent. These couplings are accounted 
Cor but, consistently with the previous discusssion, their influence is 
assumed to be constant with frequency, and their value is taken equal to 
their actual value at the roll peak frequency. With these api^roximations, 
the values of the added mass, damping, and hydrostatic terms, as well as 
their variations with speed and heading angle, are contained in Appendix A. 

Figures 3.9, 3.10, and 3.11 show the overall result of the approxi- 
inations on the ship dynamics. The highly-tuned second-order behavior of the 
roll transfer function is obvious, and the approximations appear to be 
very good. The higher frequency dynamics are neglected, but the approxi- 
mations are in a sense conservative. 

The derivation of a state-space model from the rational transfer 
functions obtained for the sea and the forces is straightforward. For 
the ship dynamics, some care has to be taken. 

Taking the matrices A^ and constant, the differential equation: 




(3.13) 


is of order 6. By selecting R = (M+A^^) P = (M+Aj^) and 

Q (MfAj^)~^Cj^, the following state-space representation can be obtained: 


r • 1 


f 

- 


■ • 


* • 

^h 


-p 

-Q 



+ 

R 



I 

0 




0 

• 


(3.14) 


The state-space representation includes 6 states (3 motions and their 
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ROLL TRRNSFLR FUNCTION 



OJg (RAD/s) 

Figure 3.10 ; Roll Transfer Function and its Approximation 
(U=l?.5£t/s, (|)=45deg) 



YRW TRRNSFKR FUNCTION 



COg (rad/s) 

Figure 3.11 ; Yaw Transfer Function and its Approximation 
(U=15.5ft/s, f-=45deg) 
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derivatives) . 

When combininq this model with the model for the forces, two pole-zero 
tMnc<5i Jutions occur. The sliip dynamics include two [hjIos at. the oriqiu, 
wtiich are due to the zero spring constant for sway and yaw, while the forces 
dynamics have each two zeros at the origin, as previously mentionncid. These 
pole-zero cancellations lead to a non-mlnimal order system, and to instabi- 
lities in numerical simulations. By simple algebraic manipulations, detailed 
in Appendix A, this problem can be resolved, and a 4th order model of the 
shiiJ dynamics is obtained, including the sway, roll, yaw motions, and only 
ttie roll derivative. 


.? . 5 Overall Ship Model 

The overall ship model is written in state-space form ; 








The state vector contains 16 states ; 


^1 •• 

wave elevation (ft) 


"^3 = 

wave elevation der l,vative 

(ft/s) 

0 

X 

X. : states related to the 
6 

sea dynamics 


/ sway force (tons.s) 


Xy : 

sway force (tons) 



/ roll moment (tons.ft.s) 


^10 •• 

roll moment (tons. ft) 


^11 ' 

/ yaw moment (tons.ft.s) 



(3.15) 
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moment (tons, ft) 

: sway displacement (ft) 

Xj^^ : roll rate (rad/r.) 

Xj^g ; roll angle (rad) 

Xj^^ : yaw angle (rad) 

The structure of the matrix A ,efined on Table 3.1, and the 

numerical values for the condition : H=10ft, 0)^=>0.72rad/s, U=i5.5ft/s, 
and <j)=45degrees are given in Table 3.2 and Table 3.3 (for any condition, 
see Appendix A) . 

The vector ^ is a white noise vector whose only non-zero element 
is the 6th row element. Its spectral intensity is 

3 . 6 Additional Comments 

Weather conditions of zero visibility and sea state 5 are considered 
in this study. Sea state 5 corresponds to a significant wave height of 10ft 
and a modal frequency of 0.72rad/s, for fully developed seas. The ship 
speed will be assumed to be 15.5ft/s (about lOkts) and the wave heading 
angle 45 degrees, A second sea condition will be sometimes used. It 
corresponds to the condition numbered condition 4 in 1 11] , with significant 
wave iieight 12ft and modal frequency 0. B07rad/s (this is the less common 
case of decaying seas). 

The ship model poles and the rms motions values (at the reference 
point) are given in Table 3.4. The model poles are also shown on Fig. 3. 12, 
for the case H=l0ft, ih^=0.72rad/s. Six poles correspond to the sea spectrum. 
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Six other poles correspond to the force dynamics (2 for each motion). The 
four remaining poles correspond to the ship dynamics : there are in fact 
six poles (two for each motion), but two pure integrators in sway and yaw 
were eliminated with zeros at the origin. There is H Very lightly damped 
mode with frequency close to 0.5 rad/s, corresponding to the roll oscilla- 
tion mode. Another lightly damped mode of low frequency appears, but its 
contribution is small, due to the presence of zeros at the origin. 

For the VTOL landing problem, the motions of interest are the motions 
at the landinc, pad. The sway motion at the landing area is composed of the 
sway at the center of the axes, plus contributions from the roll and the 
yaw angular motions. In the ship coordinates, the landing pad is located 
at ; 


^SLP (3.16) 

z 34 ft (3.17) 


The sign conventions for the aircraft are also different than those for the 
ship, so that the output matrix providing the sway, roll and yaw motions of 
the landing pad in aircraft coordinates is : 


'SLP 


000000000000 -1 0 34 127 

0000000000000010 
0 00000000000000 -1 


(3.18) 
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) . 7 Sunanary 

In this chapter, a ship model was derived from hydrodynamic data 
obtained from the geometry of the DD963 ship hull. The ship motion riiodel 
is divided in a sea model, and a ship model. The sea spectrum appears to 
be concentrated in a narrow frequency band centered around 1 rad/s, so 
that all motions have relatively limited power spectra. Roll motion 
especially appears to be highly concentrated in frequency, as it behaves 
like a lightly damped second-order system. 

Some approximations had to be made in order to obtain a finite 
dimensional model of reasonable order. It was demonstrated that these 
approximations were very good in the frequency range of interest. 

The resulting model, expressed in state-space form, and its salient 
characteristics were discussed. 
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CHAPTER 4 


SHIP MOTION ESTIMATION AND PREDICTION 


4 . 1, Ship Kalman Filter 

Tho Kalman filter role in an LQG design was indicated in chapter 2. 
Here the Kalman filter tasks are the following : 

- reconstruct the states that are not directly measurable, as for 
example the sea states (this is essential to the prediction of the ship 
motions) 

“ provide optimal estimates of the states, including those that are 
measurable, but are affected by noise. In this sense, the estimator is 
really a filter that filters the noise affecting the measurements. 

The measurements of ship motions are affected by noise that is not 
simply the instrument noise, but is also caused by the structural vibrcT- 
tions. These vibrations can be quite significant, and we want these hig!i 
frequency motions to be filtered by the Kalman filter, and to keep the 
filter poles within the range of frequencies which is significant of 
ship motions. 

The meaS’,. V '"'enl. equation is ; 




(4.1) 


It is assumed that the only measurements available are the sway displacement, 
the roll angle and the yaw angle, so that : 
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0000000 0 00001000 
000000000 0 00001 0 
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 


(4.2) 


The filter is desiyned assuming the measurement noise spectral intensity 
matrix to be : 


Q 


s 


tr: 


0.1 0 0 
0 0.0002 0 
0 0 0.0002 


(4.3) 


and the filter poles are within a radius of 1.1 rad/s. The filter poles 
and the predicted nns estimation errors are indicated in Table 4.1. 

The rms estimation errors are the square roots of the diagonal 
elements of the covariance matrix of the error : 


ee 


E(^^.e^) 


Having : 


-S S-S ^ 








= AgXg + (Xg 




The error is defined by 


(4.4) 


(4.5) 

(4.6) 

(4.7) 



(4.8) 
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Table 4.1; Ship Kalman Filter Poles and ,Rms Estimation Errors 
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The error is governed by the differential equation : 


^ ■ <*s-Vs>Ss - k - 


(4.9) 


SO that is the solution of the 16th order Lyapunov equation : 


>''s-"s=s>'’,e * ^ »sSs"I - “ 


(4.10) 


figures 4.1, 4.2, and 4.3 illustrate time simulations of the ship 
motions and of the estimated motions by the Kalman filter. 

4 . 2 Sensitivity of the Estimation Error to Parameter Uncertainty 

A study of the sensitivity of the estimation error to parameter 
uncertainty indicates those model parameters that are important in the 
estimation process. 

The error equation of the Kalman filter with incorrect model is given 
in [ 16] . Assuming the correct model to be : 






and the filter designed with A* and C*, so that 


4 -- A|»s 4 


The error equation is ; 


^ - <A*-HgO*)eg f (A5-Ag-Hg(C*^Cg))Xg - " »S^ 


(4.11) 

(4.12) 


( 4 . 13 ) 


(4.14) 
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Figure 4.2; Roll Actual and Estimated Motion 
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DonoUinfj : 
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P 
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“ B(eg.i3g) 




“ E(Xg.^) 


(4.15) 

(4.16) 

(4.17) 


and : 


P' - 


p 

P^ 

ee 

xe 

P 

P 

xe 

XX 


(4.18) 


''s 


0 A„ 


H ' 

“s 







(4.19) 


(4.20) 


The matrix P' is the solution of the 44th order Lyapunov equation : 


A!P' + P'A’’^ + H* =0 (4.21) 

b S b 

This provides the error covariance matrix and the actual rms estimation 
errors. Table 4.2 shows the result of the sensitivity study. The most 
important parameter is the sea modal freq\iency and this indicates the 
importance of an accurate estimation of this parameter aboard the ship, 
in a real-time application. 

The influence of systematic measurement errors is studied by using 
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HEA ll-=10ft 0 ) =0.72rad/s 

m 



SWAY (ft) 

ROLL ( ° ) 

YAW ( ° ) 

Kins motions 

0.G12 

4.56 

0.227 

Ni-'minal lorors 

0.241 

0.56 

0.0776 

Pcirameter cliiuuied : 




U--20ft/s 

0.245 

0.568 

0.0963 

u =-0.5 2 rad/s 
m 

0.314 

0.60 

0.0858 


0.296 

0.624 

0.112 

C (sway)-=0.9 

0.255 

0.586 

0.081 

C (roll) .=0.9 

0.247 

0.708 

0.0808 

C (yaw) '-0.9 
s 

0.242 

0.56 

0.0777 

( sway) =0 . 

0.518 

1.21 

0.1408 

C (roll)=0. 
s 

0.376 

4.08 

0.158 

C_^(yaw) =0. 

0.242 

0.563 

0.0785 


Table 4.2 ; Sensitivity of the Estimation Errors to Parameter Uncertainty 
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a calibration factor in the matrix. In the case of a 10% error, the 
most significant influence is obtained for a roll measurement error. 
Clearly, the Kalman filter uses mostly the roll measurement in the 
estimation and improves its estimation of the other motions through the 
known couplings between the motions. In this regard, the yaw measurement, 
which for this application is very noisy, has little influence on the 
filter performance. In the case of a calibration factor 0 (indicating 
a disconnected measurement), significant errors result, especially for 
roll . 

4 . 3 Ship Motion Prediction 
4.3.1 ^o^roduc^io^r^ 

Real-time prediction of ship motions is currently a subject of great 
interest, with applications not only to aircraft landings (or helicopter, 
or VTOL landings), but also to many other operations such as ship motion 
compensation and cargo transfer. Recently, a method has been proposed, 
that uses time series analysis [ 17] . Such method does not require any 
precise ship modelling : it considers the ship as a black box. 

In this research, we chose to take advantage of the available know- 
ledge on ship motions, and to derive a better model of the ship a priori, 
that can be used for estimation and prediction of the ship motions. 

Note that although the prediction of ship motions is of primary importance, 
to any piloted landing application, it is useless to an automatic landing 
using the LQ/LQG methodology, because the knowledge of the predictable part 
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of future ship motions is in fact compietely included in the ship states 
and is fedforward by the LQ/LQG controller. 


4.3.2 jPi^ed_icjtion_w^th Cor^r£c^ 
Given the system ; 


4 = ^ 

The optimal predictor is simply given by 


The error equation is : 


(4.22) 


( 4 . 23 ) 


2s " - % 


% - AsSs - £s 


(4.24) 

( 4 . 25 ) 


The cc v^ariance matrix of the error is time-varying, and is the solution 
of the linear differential matrix equation : 


P 

ee 


T 

A„P + P 
S ee ee S 




( 4 . 26 ) 


In a first approach, we assume that the state at t=0 is perfectly known, 
i. e. : 


^(0) = 2^(0) (4.27) 

For t-x», since A is stable, the estimate goes to 0. This reflects the 
fact that the ship motion is stochastic and that the knowledge of the state 
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at any time does not imply any knowledge about the state at a distant 
time in the future. Consequently, the error eventually reaches the rms 
value of the mo cion. 

The evolution of the rms error, normalized by the rms motion is 
shown in Fig. 4.4. Roll appears to be the easiest motion to predict. 

Under perfect conditions, good prediction could be obtained for up to 
10 seconds. This reflects the highly-tuned shape of the roll power 
spectri.im, which is mostly a second-order very lightly damped oscillator. 

The above is a significant result for landing purposes, because roll 
motion is really one of the most important motions to predict. From the 
lateral motions, the yaw motion is very small, while the sway motion 
consists mostly, at the landing pad, of a roll induced motion. 

Sway and yaw are much harder to predict, as shown in Fig. 4. 4, which 
reflects the much broader power spectrum of these motions. 

The oscillations in the prediction covariance can be explained by a 
simple derivation for a second-order system. 

Assuming a second-order system : 


^ - A}^ -f £ 


where : 



(4.28) 


(4.29) 


(4.30) 
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Figure 4 . 4 : Rms Prediction Errors (H=i0ft, u) -0.72rad/s) 


15.00 



Thu covariance matrix px'OpcUjation equation is : 


• 111 
P «= A P + P A + H 


(4.31) 


with : 


'll 

^12 

'l2 

^’22 


The equation (4.31) can be rewritten : 


‘’ll 


0 

2 

0 


‘^ll 


0 

‘’l2 


-b 

-a 

1 

• 

‘’12 

+ 

0 



0 

_ 

-2b 

-2a 


"22 

IB 


1 


(4.32) 


(4.33) 


This is a linear differential matrix equation. The eigenvalues of the 
above 3x3 matrix are : 


A = -a 


(4.34) 


A -ai /a^- - 4b 


(4.35) 


Thu eigenvalues of the original system are : 

A = (-a i /a^“~4b)/2 (4.36) 

In other words, the linear system modelling the covariance matrix 
propagation has eigenvalues which are equal to twice the original system 
eigenvalues. This explains the oscillations of the prediction errors which 
are at a frequency double that of the ship motion peak frequency. 
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Time simulations of the ship motions and of the prediction estimates 
are reproduced in Pigs, 4.5, 4.6, and 4.7. The prediction starts at t=40sec. 
The previous remarks on the predictability of the motions are easily so?en 
in these plots. 

In an actual realization, the ship states are not directly available 
and must be estimated. Pig. 4.8 shows the structure of the estimator/pre- 
dictor. Prediction can be seen as an extension of filtering, with the 
measurement branch broken. In an actual application, the Kalman filter 
would be running in real-time, continuously, while the predictor would 
compute the predicted motions for a specific time At ahead, in parallel. 
This is indicated schematically in Fig. 4.9. 

Some simulations are shown in Figs 4.10, 4.11, 4,12 and 4.13. The 
Kalman filter provides estimates at t=40s that are used to predict the 

ship motion. Good prediction is still o; i led for roll, while yaw becomes 

almost unpredictable. 

4.3.3 

The actual implementation of the estimator/predictor described above 

will require the use of a digital computer and, hence, the discretization 

of the differential equation. Although the Kalman filter requires a large 

number of computations at relatively small time steps (with measurement 

update) , the particular structure of the predictor makes the additional 

computational load very small. Given the optimal estimate of the Kalman 

filter at time t, x^(t) , the optimal estimate of the predictor at 
—fa 
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Figure 4.7 ; Vaw Actual and Predicted Motion 
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Figure 4.8 ; Estimator/Predictor Structure 
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Figure 4.9 : Estimator/Predictor Practical Structure 
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Fi gure 4.11 ; Rol.l Prediction with Noisy Kaasurements 
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Figure 4.1 2: Yaw Prediction with Noisy Measurements 
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t. i.me t,+Ab is simply : 

A,, At. 

Xg(tiAt) -• o ^(t) (4.37) 

A^At 

Having computod e , and (;m\sidoriru) that we are only interested in 
predicting the motions (and not the other states), this operation results 
in multiplying a matrix (3x16) by a vector (16x1). These are the only 
computations required if we are only interested in predicting the motions 
Al seconds ahead (and not their evolution in between) . In this case* there 
is no need tc integrate the differential equation, as for the Kalman 
filter. 

In the case when the complete time evolution is desired, and that the 
differential equation has to be integrated numerically, the same equation 
(4.37) provides the discreti'^ed equation corresponding to the differential 
equation : 






(4.38) 


It is very tempting to approximate this equation by the simple diffox'ence 
equation : 


S (ttAt) = ^(t) + At(A ^(t)) 


which is thi. scune as approximating 
A At 

e I 4- A„At 

S 


(4.39) 


(4.40) 


for cl small At. 
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However, due to the presence of the very lightly damped roll mode, large 

errors can occur even for sanall At, compared to the time constants of the 

system. In fact, if At is not very small, the mode becomes even less damped, 

AgAt 

as indicated in Appendix B. More precise approximations to e are 
advisable for the implementation of the predictor (see for example [18] ). 
Such problems do not appear for the Kalman filter which does not have this 
very lightly damped mode. 


4.3.4 Pr^<iic^ion_w^th Wrong^ 

In a previous subsection, the important inf luence of the sea modal 
frequency on the performance of the Kalman filter was assessed. This para- 
graph addresses the same question for the predictor. The calculation of 
the errors is a little more complex in this case. 

Given the system ; 


is “ 


and the predictor with wrong model 


4 ■ *S is 


The error is given by 


(4.41) 


(4.42) 


% - l''|-''s>is + ''SSs ' 5s 


(4.43) 
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The global system is : 



(4.44) 


(4.45) 


(4.46) 



P' is the solution of the linear differential equation 


* T 

P" = A' P' + P' A' + S' 
S S S 



(4.47) 

(4.48) 

(4.49) 

(4.50) 

(4.51) 


The original system (the ship) is assumed to be in steady-state, so that 
P *0, and P is the solution of the Lyapunov equation ; 

XX XX 
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0 


(4.52) 


A^P + P A^ 
S XX XX S 


S 


Replacing P by its steady-state value (given by (4.52)) in equations 

XX 

(4.50) and (4.51), we obtain a 32nd order linear differential matrix 

equation in P and P , with initial conditions P (0) = P (0) = 0 
^ xe ee xe ee 

(the initial state is assumed to be perfectly known) . 

Figures 4.14, 4.15 and 4.16 show the results of the integration of 
this equation, in the case of an error in the sea modal frequency. The 
filter was designed with o)^=>0.72 rad/s, while the actual value was 
assumed to be 0.52 rad/s. The degradation in performance is again signi- 
ficant and the roll prediction time is about divided by two. This 
indicates (again) the importance of the estimation of the sea modal 
frequency in real-time applications. 

4.4 Summary 

In this chapter, we addressed the important questions of ship motion 
estimation and prediction. First, a Kalman filter was designed, whose task 
is to filter the noise in the measurements (mostly due to the ship structur- 
al noise), and to provide estimates of the states that are not available for 
measurement. These estimates can then be used for optimal prediction of the 
ship motion, and in the feedforward path of the LQG controller for the VTOl 
landing. 

A sensitivity analysis showed the relative importance of the model 
parameters, and the sea modal frequency appeared to be a significant 
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Figure 4.14; Sway Rms Prediction Error with Wrong Model 
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parameter. 


The prediction of ship motions, assuming perfect knowledge of the 
initial state, was studied, so that lower bounds on the prediction errors 
were obtained. Due to its concentrated power spectrum, the roll motion 
turned out to be the easiest to predict (optimally 5 to 10 sec. prediction 
time). The sea modal frequency was shown to be an important parameter for 
all motions, which indicated the need for an accurate esti nation of this 
parameter in a real-time application. 
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CHAPTER 5 


AIRCRAFT MODEL 


5.1 Introduction 

Since the beginning of VTOL aircraft technology development, a large 
number of configurations were conceived, and sometimes built and tested. 

The most studied and extensively tested VTOL is probably the AV-8A Harrier. 
Another VTOL aircraft which has been extensively studied is the Lift/Cruise 
Fan V/STOL Research Technology Aircraft (RTA) . For this aircraft, complete 
simulator programming data is available in ( 19l . In the former study of 
the longitudinal controls required for VTOL landings on destroyers 111 , 
a linearized model, written in state-space form, was derived from the data. 
It includes both longitudinal and lateral motions. In addition to the rigid 
body equations of motion and the contributions of the fan forces and 
moments, this model also accounts for ram drag forces and moments, and 
internal momentum effects (gyroscopic terms) due to the rotating engines 
and fans. Neglected are the aerodynamic effects, the ship airwake turbulence, 
and the ground effects. The actuators are modelled by first order dynamics, 
whose time constants are to be selected by the control system designer. 

5. 2 Lateral Motions Model 

To a first approximation, the longitudinal and lateral motions of the 
aircraft can be decoupled. This corresponds to neglecting the cross-coupling 
terms between roll, yaw, and pitch due to gyroscopic effects from the 
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engines and Cans. The effect of this approximation on the overall system 
performancQ and stability margins is an interesting issue which is not 
ac' dressed in this thesis, but is probably worth a subsequent study. 

The general ship-aircraft configuration is shown in Fig. 5.1. The 
lateral motions, and the sign conventions used, are indicated: they are 
sway (lateral translation) , roll ana yaw (angular motions) . The 

definition of the controls requires some care. Strictly speaking, there 
are 9 variables on v/hich the controller can act: the 3 values of the 
thrusts at the* fans, the 3 values of the longitudinal deflections of 
these thrusts, and the 3 values of the lateral deflections. However, 
these controls are not independent: for example, deflecting the thrusts 
T^ and T^ lat<?/,ally and in opposite directions would result in counter- 
acting forces, and, consequently, in an inefficient use of the controls 
available. Considering reasonable use of the controls, we find three 
independent controls for the ]a teral motions (corresponding to the three 
degrees of freedom) ; 

1) an equal lateral deflection of the aft louvers, denoted by 2 

2) an exchange of thrust from T^ to T^, denoted by 6 t^^ 2 (this 
can easily be done, as the engines 1 and 2 are mechanically coupled) 

3) a lateral deflection of the front louver, denoted by 6a^. 

These control variables are grouped in a vector denoted by £. 

In the original work by McMuldroch [ 1] , the vector inp'jt was chosen 
in a slightly different manner. It is denoted by ij in this thesis. 
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Its components are: 


1) 6y : an equal lateral deflection of a?l the louvers, to produce 

c 

(mctly) a sway acceleration 

2) 6^^ I an exchange of thrust from to to produce a roll 
^‘cceleration 

3) : an opposite lateral deflection of the front and aft louvers, 
o 

to produce a yaw acceleration. 

In this thesis, we keep the McMuldroch formulation, but the control 
weights in the quadratic cost and the robustness criterions will be 
expressed in terms of the vector £ , which is the actual physical control 
input. The transformation from £ to u and u to £ is straightforward. 
Figure 5.2 indicates the structure of the aircraft model. 





Figure 5.2; Aircraft Model Structure 


C - 


- 96 - 



The numecical values of tlie matrices are given in Table 5.1. The units of 
the model are: 

- ft for y 

- ft/s for y. 

A 

- rad for ®“l,2' *^^3 ' *^^ 0 ' '^'^c 

- rad/s for ii»^ 

■ fraction of the nominal tlirust for 6Tj^ 2 nominal thrust of ^ 

is 9057.5 lb) and for 6(|i 

c 

'Die only difference between the model Indicated here and the model 
given in [1], is the absence of actuators dynamics which, in[l], were 
selected as integrator dynamics. For the lateral case, the actuators 
dynamics are well beyond the significant ship motion frequencies (0.5 to 
1.5 rad/s): the thrust deflection can be done almost instantenuously, 

while the exchange of thrust which does not require any overall increase 
of power - has a time constant of about 0.1 s (a shaft-coupled configuration 
is considered here). For these reasons, and for simplicity, the actuators 
dynamics are neglected here. 

The resulting model is very simple: it contains only 6 states, 
specifically the aircraft motions and velocities. As indicated earlier, 
a simple model is desirable on the control system point of view, although 
it necessarily implies some crude approximations and, therefore, increases 
thiB need for good robustness properties. On the other hand, the interpre- 
tation of thv'' elemeivts o£ the model is easier from a simple model, as will 
be seen in the next section. 
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8LS64 

,1830 
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Table 5.1: 


Aircrait Model Values 



5*3 Model Interpretation 


The important couplings inherent to a VIOL aircraft are represented 
in the model, and it is worth taking some time analyzing them. 

A first in^wrtant coupling is a dynamic coupling between roll and 
sway. Aa the aircraft is banked, a lateral con\ponent of the thrust appears, 
which results in an important side force. This is represented by the 
A^(4,2) element (cf table 5.1). It is equal to g (32.2 ft/s^) because, 
under the small angle approximation, the side force is equal to the weight 
times the roll angle. 

The other elements coupling the velocities in the A matrix are due 
to the ram drag forces and moments, and are responsible for most of the 
aircraft dynamics in this model (especially its instability) . 

An important coupling appearing in the B matrix is the B (5,1) term. 

A A 

It indicates tliat a lateral deflection of nhe thrusts produces an important 
roll moment. Note that if this effect is not compensated for, the roll 
moment will produce a roll angle which, by the effect described above, will 
in turn produce a jway force opposite to the sway force originally produced 
by the louvers deflections. The origin of the roll moment is in the 
difference between the center of gravity and the center of thrust of the 
fans (the center of gravity is 3.12 ft higher). 

Finally, an interesting, although apparently not very significant term^ 
is the B^(6,/2) term. It tells us that an exchange of thrust in the aft fans 
produces a yaw acceleration. Although no yaw moment is produced by this 
control, a small yaw acceleration resultfii from the roll moment, due to the 
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angular difference between the principal axes and the body axes (in other 

words, due to the presence of a cross-term 1 in the inertia tensor). 

xz 

The above summarizes the important couplings present in the VTOL 
aircraft. TJiese couplings are dynamic (reflected in the matrix), as 
well as control couplings (reflected in the matrix), and are captured 
in tlie simplified, linearized model given here. 

In fact, not only is tlie system strongly coupled, but it is also 
open-loop unstable. The open-loop jxjles are shown in Fig. 5.3. Two of the 
six poles are at the origin and correspond to the pure integrators in 
sway and yaw motions. 

5.4 Summary 

In this chapter, we described the aircraft model for the lateral 
motions. We showed tlnat three independent controls can be used to track 
the lateral deck motions, but that some important couplings are present 
amongst the controls and motions (especially between sway and roll). 

The aircraft model is written in state-space form, and the state 
vector contains 6 states, namely the aircraft motions and velocities. 

T)ie important couplings of the VTOL are represented in the model, and 
it was shown to be open-loop unstable. 
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CHAPTER 6 


LINEAR QUADRATIC REGULATOR DESIGN 


6.1 Introduction 

The LQ/LQG design methodology was introduced in chapter 2. As a first 
step, it is assumed that all states of the system are available (output 
matrix C « 1), so that full state feedback can be used. Note that while 
this assumption is generally a rather unrealistic assumption, it is not 
the case for he simplified aircraft model that was obtained in chapter 5. 
Angular motions as well as angular rates are available through gyros and 
rate gyros. Similarly, estimates of position as well as velocity can be 
obtained from accelerometer measurements combined with MLS/DME measurements, 
with a high quality. In a sense, a Kalman filter is not strictly necessary 
for the aircraft. For these reasons, significant attention has been given to 
this part of the design. 

^ • 2 Choice of the Quadratic Weights 

Under the assumption of full state feedback, the only parameters left 
to the designer are the Q and R weighting matrices in the quadratic cost. 

A very natural way to select these matrices is the diagonal inverse-square 
weighting [20]. Some other methods have been proposed, as for example the 
method proposed in [21] and [22] to achieve desirable asymptotic regulator 
properties. In any case, it is interesting to note that the robustness 
properties of LQ regulators can be seriously deteriorated if a non-diagonal 
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matrix R is chosen (23] [24]. For this reason, it is wise to select a 
diagonal R matrix, or equivalently, to pick R » pi, after some scaling 
of the inputs. 

In this thesis, we decided to start with a simple diagonal inverse 
square weighting, leaving some parameters to acheive a desirable eigen- 
structure, with a careful study of their influence on the optimal root- 
locus. The quadratic cost that we want to minimize xs given by (2.25). 
As a first step, we wane to weight all the state variable, (motions and 
velocities), so that: 



The Q matrix is clxjsen diagonal: 
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where y , d) , d' are arbitrarily choeen maximal deviations of the 
'max' ^max ^max ' 

states: 


'max 


ft 


(6.5) 


<l> 


max 


'P. 


max 


10 degrees = 0.174 5 rad, 


( 6 . 6 ) 


A parameter is left in the Q matrix. For «=> 0, only the motions are 

weighted. For T 0, the motions and the velocities are weighted. They 
z 

are weighted equally if T = 1. In the VTDL landing problem, the velocity 

z 

tjracking errors can be as important as the position tracking errmrs 

(they determine the shock at the landing) , which justifies the weighting 

of both in the quadratic cost. Also, it is known that, sometimes, the 

weighting of outputs only may lead to underdamped second-order type 

systems, and that this can be improved by weighting derivatives as well 

as outputs. At this stage, the variable T is then left as a design para- 

z 

meter. 

The co,ntrol weighting R matrix is chosen diagonal at the point of 


the actual controls: 


^6a 


1,2 


max 




0 

JL 

1.2 


max 




0 

1 

i 

max 


(6.7) 
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and: 


T 

T R T 
c c c 


( 6 . 8 ) 


with 6aj^ 2 » 2 » 

' max ' riax "^max 

of the control authority. 


are arbitrarily clioeen maximal values 


fia, „ ■ 6a_ " 10 degrees ■ 0.1745 rad (6.S) 

' ^roax ■^max 


6T| „ " 0.3 (30% of the nominal thrust value) (6.10) 

X I 

max 


Another parameter of the design, p, is left here, and indicates the 
relative importance of the control versus state deviations. 


6.3 Optimal Root-Locus 

An optimal root-locus is defined as the locus of the closed-loop 
poles of the system with optiiaal LQ feedbaclc, when the parameter p in the 
R matrix is varied from “ to 0. 

It is interesting to note that this proJslem in itself does not require 
the solution of the associated Riccatti equation. The closed-loop poles are 
the left half plane eigenvalues of the Hamiltonian system [9] and[251 : 


Z n Z Z 


with: 



Z “ 


-Q 


p " 




( 6 . 11 ) 


( 6 . 12 ) 
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t*h 

This repltices the problem of solving a n order Riccatti equation (n being 
the order of the system) and a n*^^ order eigenvalue by a 2ri^^ order 

eigen faluo problem. 

The points of departure of the optimal root-locus are the stable 
open-loop poles, combined with the mirror images of the unstable oj)on-loop 
poles (symmetric position, ^bout the imaginary axis). In our case, two poles 
are at the origin, two are on the negative real axis, and two are oscilla- 
tory unstable, and are then reflected in the left"half plane. 

The points of arrival define the asymptotic behavior of the system as 
p-KO, and are somewhat harder to obtain. Some poles go to infinity along 
specific asymptotic patterns (Butterworth patterns) depending on the rank 
of Q**’’. The other poles go to the zeroes of transmission of the system 

L 

Q (sI-A) B (with the same remarks for right-half plane zeroes as for 
unstable poles when p-»o>) . The zeroes of transmission z^ are the solution 
of the generalized eigenvalue problem: 


(A-z^I) 

-o'* 



0 


(6.13) 


The system has as many inputs as outputs, so that the zeroes of transmission 
are the z.'s that make the matrix A' (defined by (6.13)) rank degenerate 
(the full rank is m+n, where m is the number of inputs, and n the number of 
states) . 

The matrix Q is diagonal, so that we can take: 
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Q 


S 


» dl«g 


1 1 1 
m«x ^max ^max 


-* max ^max ^max 


(6.14) 


and A' is given in Table 6.1. 

The matrix is three times rank degenerate at s, » -1/T . Than the row 1 

i z 

is a linear combination of the rows 7 and 10 (similarly for 2, 6, 11 and 
3. 9, 12). In other words, the inclusion of the velocities in the 
quadratic cost results in the introduction of three transmission zeroes 
in the optimal root-locuc. It turns out tliat the three remaining zeroes 
of transmission are at the infinity. The case when T »0 can be seen as a 

Z 

limiting case when the weights on the vderivatives go to zero. All the 

transmission zeroes are thei; at infinity. 

Figures 6.1, 6.2, and 6.3 show the optimal root-loci for T «0, T «»l, 

z z 

and T sQ.S and are labelled for different values of p. The structure of the 
z 

root-locus for high values of p is the same for different values of T 

z 

and is not repeated in the last two plots. 

The previous comments concerning the points of departure and cf arrival 
are easily checked. For non-zero T , t)w poles eventually reach the negative 

Z 

real axis, and three go to infinity, while the others reach the zeroes at 

-1/T . As T goes to zero, the zeroes move to infinity, and the case of zero 
z z 

velocity weights appears clearly as a limiting case from the three figures. 
Then, the asymptotic structure consists of three second-order Butterworth 
patterns. 

1 exp/^cted, the root-loci reach faster higher damping regions when 
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Table 6.1: A' Matrix Structure 





ORtGiNAL PAQS IS 
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derivative weights are included. However^ the asymptotic behavior with 
no velocity weights is a desirable one. At this point, no decision is 
made about which case is better (T »0 or T »>1). This depends on which value 
of p corresponds to a realistic design (acceptable tracking errors and 
control autiKarity) . If high values of p are reasonable, choosing T «=1 
will probably improve the design. If small values of p are reasonable, 
choosing T_ =*o may be better, as it loads to desirable closed-loop poles 
locations. 

It turns out that the ranga of p between p » 0.3 and p « 30 is a 
realistic one, as will be elaborated in the sequel. Tlie corresponding pole 
locations are satisfactory for “ 0, and this case will be considered as 
the nominal one from now on. 

Finally, an interesting characteristic of the root -loci presented 
here is their rather peculiar behavior around the 1 rad/s region. This 
appears in all three cases. For T = 0, the behavior is quite surprising, 
as one pole (one complex pair) seems to come back before going to infinity 
along the 45 degrees line. At some point, this point "slows" down as if 
it was reaching a zero. This point is found to be approximately at 
0.91 ± j 0.87 and p ~ 3. No attempt is made to justify this behavior 
mathematically, but some physical connection can be found, and will bo 
explained below. 

It can be expected that the imixjrtant couplings between sway and roll 
previously mentionned have some importance in this strange behavior. 
Actually, while a small relative change of p produces only a very small 
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itK>vement of the pole near p>« 3/ a comparatively email change of the weight 

on 6T - (6T. „ ) produces a much larger movement of this pole. It is 

' ' ^max 

clear then tliat the roll control variable is crucial in the optimal lateral 
control system design, as it is probably in any VTOL lateral control system 
design. 

6.4 Step Responses 

As a first step in the evaluation of the control system design, 
responses to deviations from equilibrium are computed and plotted. The 
aircraft is assumed to be left with a zero velocity and some position error 
at t=0. Figures 6.4, 6.5, 6.6, 6.7, show the response to a 4 ft sway 
initial error, and Figures 6.8, 6.9, 6.10, and 6.11 the responses to a 10 
degrees roll initial error. The responses to a yaw initial error are not 
significant (they require comparatively very small control authority), and 
are not reproduced here. Similarly, the deflections ^ and are 
about the same, and only 6a is slK>wn. 

X f A 

The sway error response is particularly interesting: for high values 
of p (high control cost), the roll response is important, while the fan 
deflection (which mainly produces a sway force) is small. In fact, the 
controller flies the VTOL as a helicopter: as the weight on the state is 
small, and the weight on the controls is high, the controller slowly banks 
this aircraft and uses the lateral component of the thrust to obtain a sway 
acceleration; then, after some time, it banks the aircraft in opposite 
direction to reduce the sway speed to zero, together with the sway deviation. 
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Figure 6.5 ; Roll Response to a Sway Initial Error 
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Figure 6.6: Fans Deflection Response to a Sway Initial Error 
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Figure 6.7 t Thrust Exchange Response to a Sway Initial Error 
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Figure 6.8: Sway Response to a Roll Initial Error 
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Figure 6.9; Roll Response to a Roll Initial Error 
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Figure 6.11; Thrust Exchange Response to a Roll Initial Error 
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For BRUkll p (high state deviation cost), the situation is very 
different, and a dramatic increase in fan deflection response can be 
observed on Fig. 6.6, while the roll response on Fig. 6.5 is much 
smaller. This shows that the turning point observed in the root-locus 
of Pig. 6.1 corresponds also to a ctiangc in strategy of the controller. 

A similar conclusion can be obtained from the responses to a roll 
initial error. For a high value of p (p«30), the sway induced by the 
initial roll angle is important, and it is then compensated by an impor- 
tant opposite roll (about lialf of the initial roll angle) . 

We conclude from this discussion that the quadratic optimization 
problem leads to solutions which, when understood, are very logical, 
and simply express specific characteristics of the system, which, at 
first, may be obscured by its multiple-input multiple-output structure. 
One advantage of the LQ methodology is that it leads to coupled 
controller designs that exploit the dynamic coupling phenomena. 

6 . 5 A Special Example 

To illustrate the above remarks, a limiting case is presented; it 
corresponds to p <= 10, T = 0, but with the penalty on the sway error 

Z 

multiplied by 10®. The resultant clc .d-loop matrix: 


A, CL 


\ - 


B G 
A A 


(6.15) 


and the gain matrix are given in Table £.2. 

The most interesting terms are the A, _ (4,2) and the G.(2,l) terms. 

n t Vtflj A 
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0. 

0. 

0. 

1. 

0. 

0. 

0. 

0. 

0. 

0. 

1. 

0. 

0. 

0. 

0. 

0. 

0. 

1. 

-32.96 

-0 . 11 

0.0155 

-25.55 

-0.3263 

0.0021 

44.77 

-12.04 

-0.4046 

2.634 

-5.656 

-0.4832 

17.93 

-0.240B 

-1.175 

1.279 

-0.3219 

-1.539 


- 

10.2360 

1.0034 

-0.00005 

0.7914 

0.0009 

■‘0.0012 

-0.0637 

3.5568 

0.0071 

0.1737 

1.2153 

0.0065 

-1.9995 

-0.0093 

0.1443 

0.1462 

0.0009 

0.1711 

Table 6.2: 

Closed-Loop Matrix and 

Feedback 

Gains for 

a High 


Sway Error 

Penalty 






Figure 6.12 ; Exact Compensation of the Lateral Force Due to 
a Roll Angle with a Thrust Deflection 
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Compurod to the original A (4,2) J2.2, the term A, (4 , 2) »-0.11 is 

” iK f Oli 

reduced almost to zero; we recall that this term expressed the lateral 
force due to a roll angle of the aircraft. The weight on a sway error 
being very large, the controller compensates this effect very logically 
with a deflection of the thrusts in opposite direction than the roll 
angle. This is accomplished by tlie G^(l,2) term, which, in the limit, 
would be equal to 1 (see Fig. 6.12). 

^ ^ Ship Mot io n Tracking 

Up to now, the VTOL controller was designed as a stabilizer, or 

zero state tracking controller. As indicated earlier, the feedback 

structure of tJjo controller is independent of the specific reference 

signal to follow, i.e. the ship motion at the landing pad. The ship 

itkJdel is required in the definition of the forward gain matrix value (GJ. 

s 

The importance of the ship model is mostly in the computation of the 
rms tracking errors and controls, which define the performance of the 
controller. 

,tn chapter 2, we denoted by w and w_ the vectors of the aircraft 

A — b 

and ship variables that we want to track. They are the aircraft and ship 
motions (at the landing pad) , so that; 



with 


W- = 

A 


1 0 
0 1 
0 0 


0 0 
0 0 
1 0 


0 0 
0 0 
0 0 


( 6 . 18 ) 


W 


s 


c 

SLP 


( 6 . 19 ) 


The tracking errors are denoted w. c' that 


5^A-S " 




W- _ X 
A-S — 


( 6 . 20 ) 


where : 


W = ( -W 
'A-S ' S 


W. ) 


( 6 . 21 ) 


X 



( 6 . 22 ) 


The overall system equations are: 


X = Ax + Bu + £ 


( 6 . 23 ) 


c = T u 
— c — 


where: 


A 



0 


0 



( 6 . 24 ) 


( 6 . 25 ) 
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A 

5 * 


(6.26) 


(6.27) 


The spectral intensity of ^ is: 




The LQ controller is designed, using a Q matrix: 


Q = ^ Q w, ^ 

^ A-S A-S 


T T 

Q -W„ Q W, 

S *w S S A 

T T 

~K Q w,, w;: q w„ 

A w S A w S 


J 


(6.28) 


(6.29) 


where Q is the same matrix as the one used in 6.2 and defined by (6.4). 

■ w 

Similarly, the matrix R is chosen as defined by equations (6.7) and (6.8). 
The control law is : 


u = -Gx >= 


G Jc 


(6.30) 


In this state-space framework, the rms values of the tracking errors 

and of the controls c_ are very easily obtained. The covariance of the sfates 
deviations is obtained by solving the Lyapunov equation: 

(A-BG)"^? + P(A-BG) +5 = 0 (6,31) 
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where; 


P = E(x X*'*') 


(6.32) 


The rms tracking errors and rros controls are the square roots of the 
diagonal elements of the covariance matrices: 


sl-s> = Vs "I-s 


T T T 

E(c C ) = T G P G T 
c c 


(6.33) 

(6.34) 


Tables 6.3 and 6.4 contain the results obtained for the two sea conditions 
discussed in chapter 3. Table 6.5 contains the closed-loop poles locations 
corresponding to the differ-»nt values of p. 


SEA H=12ft OJ =0.4807 rad/s 




Sway 

(ft) 

Roll 

(deg) 

Yaw 

(dog) 

Controls 

Ship motion 
at landing pad 

7.155 

12.64 

0.373 

®“l,2 

(deg) 

«h,2 

(%) 

«“3 

(deg) 

Tracking 

errors 

A/C driving 
noise incl. 

P=30 

1.124 

14.98 

0.414 

0.26 

2.04 

0.28 

P=3 

0.971 

13.35 

0.203 

1.93 

3.06 

m 



p=0.3 

0.465 

6.36 

0.099 

9.15 

15.35 

9.16 

P-0.3 

0.801 

6.41 

0.329 

9.72 

16.74 

9.37 


Table 6.3; LQ Controller Performance (decaying sea) 
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SKA ll-lOft U) ^0,72 rad/s 

m 


Sway 

(ft) 

Rol 1 
(deq) 

Yaw 

(deq) 

Coni I'-ols 


Ship motion 
at landing pad 

2.i)hi 

4.556 

0,227 

*“l,2 

(dog) 

®‘i,2 

(M 

6a3 

(deg) 

Tracking 

errors 

A/C driving 
noise incl. 

p=30 

0.517 

5,405 

0. 172 

0,10 

0,12 

0,10 

p=3 

0.444 

4.830 

0.082 

0,70 

1,10 

0,71 

P==0. 3 

0.227 

2.322 

0,044 

3,34 

5.61 

3,34 

P--0.3 

0.691 

2.480 

0, 317 

4.72 

8,78 

3,92 


Table 6.4 ; LQ Controller Performance (fully developed sea) 


Cl osud- 

... 1 ’• — 

loop polos locations 

p-30 

-0.598 i j 0.585 -1.038 i j 0.611 -0.642 t j 1.054 

11 

Q. 

-0.906 ± j 0.87 -1.054 ± j 1.039 -1.528 1 j 1.552 

p=0,3 

-1.064 t j 1.036 -1.869 ± j 1.863 -2.792 ± j 2.803 

i 


Table 6.5 ; LQ Controller Closed-Loop Poles Locations 
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For p B 0.3, good tracking is obtained, with reasonable control 
authority (for H «12 ft, = 0.4807 rad/s, some higher value of p may 
however be appropriate) . The value of p » 0.3 was chosen as a nominal 
value for the control system. The closed-loop poles are also at desirable 
locations for a VTOL control system. 

For this nominal design, the effect of wind turbulence is also 
indicated. The results are obtained by introducing an aircraft driving 
noise matrix 5^, and are also indicated in Tables 6.3 and 6.4. 

Values for wind disturbances are found in the form of Dryden spectra 
in [19]. The derivation of the matrix from this data is included in 
Appendix C. 

As for the aircraft controller step responses, the values of the rms 

tracking errors and controls give us some interesting indications about the 
physical aspects of the problem. The yaw tracking errors are very small, 
compared to the maximal values given in the quadratic cost. Similarly, 
the rms deflections 6 a, „ and 6 a are very close. This indicates that 
yaw tracking is not at all a problem for the aircraft. The control 
authority is clearly more than sufficient to track the ship motion 
(which is very small anyway), or to reject wind disturbances. 

On the other hand, roll tracking errors are much higher than the 
sway and yaw errors (compared proport ionnally to the maximal values given 
in the quadratic cost) . For p = 30, the rms error is even larger than 
the rms ship motion. This probably indicates that roll is the "least 
easily controllable" state of the VTOL aircraft, but it also reflects two 


- 129 - 


basic contradictions faced by the controller in the tracking of the sway 
and roll motions. The first comes from the ship motion at the landing pad. 
From the sway at the landing pad, a large part is coming from the roll 
motion, due to the difference in altitude between the landing pad and the 
center of rotation of the ship. It has been seen earlier that an easy way 
for the aircraft to track a sway reference is to roll the aircraft as 
a helicopter would do. However, the roll angle required to follow the 
ship sway motion by this mean is precisely opposite to the ship roll angle 
that has produced the sway motion of the landing pad. This is illustrated 
schematically on Fig. (;i.l3. 

The second adverse effect has been mentioned previously and is 
illustrated on Fig. 6.14. It is shown that a roll moment is induced by 
a lateral deflection of the thrusts (term B (5,1)). It is opposite to the 

g\ 

ship roll motion. 

It is possible to improve the roll tracking by increasing its penalty 
in the Q matrix. Considering the physical problems mentioned here above, 
this will probably result in a large increase in control authority. 

Moreover, precise roll tracking is not necessary, and not even desirable. 
For the same reasons as the ones mentioned here above, precise roll 
tracking would result in large lateral accelerations at the pilot location 
and this would probably be unacceptable. 

Finally, it can be noted that only the influence of wind gusts was 
considered here, while the mean wind was neglected (together with the 
aerodynamic effects, as in f 1] ) . The influence of the mean wind will be 
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Figure 6.13 ; Adverse Sway/Roll Coupling in the Ship Motion 
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to produce a steady-state tracking error whose amplitude depends on the 
amplitude of the mean wind and on the feedback gains. This can easily be 
compensated for, if necessary, by introducing some integral compensation 
(for example, by augmenting the system and recomputing the feedback gains). 

6.7 Time Domain Simulations 

Some time domain simulations are illustrated in Figures 6.15 to 6.27. 
Pigs 6.15 to 6.19 correspond to the case p ■ 0.3, witlwut aircraft driving 
noise (wind gusts). As expected roll tracking is the worst of all, and the 
aircraft roll motion is systematically smaller than the ship motion. 

Figures 6.20 to 6.22 illustrate the tracking in the presence of wind 
disturbances. As expected, the yaw tracking is the most severely affected, 
but the aircraft motion remains small . 

Figures 6.23 to 6.27 show the case p = 30, without wind disturbances. 
While sway tracking remains good, the roll tracking is simply out of phase 
most of the time. This demonstrates again the "helicopter" behavior of the 
controlled aircraft for high values of p. 

6.a> Summary 

In this chapter, we described the design of a linear quadratic (LQ) 
regulator. This design assumes perfect knowledge of the states at any 
moment. For the aircraft part, the assumption is not unrealistic (espe- 
cially if a good navigation system is used aboard the aircraft) , as all 
the str'es are available for measurement, and the only role of an optimal 
estimator is to filter the noise, and provide optimal estimates on the 
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6.15 ; Sway Tracking (p=0.3) 
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Figure 6.17t Yaw Tracking (p=0.3) 
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Figure 6.2l !Ro3,l Tracking with Wind Disturbances (p«0.3) 
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Figure 6.22; Yaw Tracking with Wind Disturbances 
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basis of the knovm dynamics of the aircraft. 

For the ship however, a Kalman filter is necessary to reconstruct 
the states that are not available, and the results indicating the rms 
tracking errors for a given control authority must be interpreted as 
lower bounds. 

Due to the decoupling of the ship and of the aircraft in the LQG 
design, the optimal closed-loop poles locus could be plotted independently 
from the ship motion. It was shown that the root-locus can be shaped by 
the introduction of velocity error penalties in the quadratic cost, and 
that some peculiar behavior of the root-locus could be explained on the 
basis of physical arguments. 

The responses of the controlled aircraft to deviations from equili- 
brium were plotted, and the important couplings between the sway and roll 
were shown to have a significant influence on the behavior of the 
controlled aircraft. For high values of the control penalty, the controller 
appeared to fly the VTOL like a he.licopter, while, for smaller values, it 
used much more the thrust deflection capability to provide lateral 
accelerations . 

The difference in strategy for high and small values of the control 
penalty appeared to have a significant influence on the errors achieved in 
the tracking of the ship landing pad motion. It was shown that the 
controller faced tiao important contradictory pheneroena in the tracking of 
the roll motion, especially due to the sway motion induced at the landing 
pad by the roll motion. Although precise roll tracking is not necessary. 


the high values that the roll motion can reach in high sea states make 
imperative the limitation of the tracking errors (to avoid a contact of 
the wing tip at the landing for example), and the constraints present in 
the simultaneous tracking of the sway and roll motion are basic 
limitations in the tracking of the lateral motions. 


CHAPTER 7 


LINEAR QUADRATIC GAUSSIAN REGULATOR DESICM 


7.1 Introduction 

Having dasigned the linear quadratic regulator, the next step to the 
final control system is the optimal estimation of the system states. As 
all the states are not available for feedback, and as those which are 
measurable are affected by measurement noise, a state estimator is necessary 
and, in the general assumptions of the LQG methodology, this state estimator 
is a Kalman filter. The overall control system is then composed of a Kalman 
filter whose state estimates are multiplied by constant gains, determined 
as in the LQ problem. 

It has been previously indicated tliat the estimation problem can be 
reasonably decoupled. The ship motion estimation problem has been addressed 
in chapter 4, so that this chapter will concentrate on the aircraft motion 
estimation problem, and on the overall system performance. 

7.2 Navigation Systems 

The accuracy of the navigation systems used aboard the ship and aboard 
the aircraft will be a determining factor in the performance obtained in the 
tracking of the ship motions. 

Some previous studies [26], [27]^ have addressed the navigation problem 
of VIOL landings on destroyers. An inertial measurement unit aboard the 
ship provides measurements of the ship motions that are data-linked to the 
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aircraft (In our configuration, the full ship state, estimated by the 
ship Kalman filter described in chapter 4, is data-linked to the aircraft). 
Aboard the ai?:craft, the combined use of an inertial measurement unit, 
microwave (scanning beam) landing system, and distance measuring equipment 
(MLS/DME), leads to estimates of the aircraft position and attitude. The 
aircraft navigation system considered in [26], [27] is a low to medium 
quality navigation system, so that even if the whole aircraft state is 
available for measurement, some improvement can be gained by an aircraft 
Kalman filter that estimates the aircraft total state using the model for 
the aircraft dynamics and control effectiveness. This is the subject of 
section 7.3. 

However, it should be noted that with the use of a good quality 
aircraft inertial measurement unit, much laetter estimation accuracy can 
be achieved. The appropriate Kalman filter formulation in this case uses 
an error state formulation, where the state variables are the very low 
frequency errors in the indicated position , velocity, and attitude of the 
i.tertial system. The radio measurements (MLS/DME) are used to estimate 
these low frequency errors. The precise structure of such a navigation 
system is not studied in this thesis, which concentrates on the control 
aspects of the VTOL landing problem. If the accuracy and response time of 
the navigation system is such that the aircraft navigation errors and the 
time lags are negligible, an aircraft Kalman filter, using the known 
dynamics of the aircraft, is not necessary, and the control loop around 
the aircraft is equivalent to the LQ control Soop. The robustness of the 
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control system is then the same as the robustness of the LQ system (cf 
chapter 8), and the degradation in performemce, as compared to the LQ 
controller performance indicated in chapter 6, is due to the ship Kalman 
filter, and the estimation of the ship states that are not available for 
measurement. The reduction in performance is small, however, as will be 
indicated in section 7.4, 


7-3 Aircraft Kalman Filter 

First, we assume that the aircraft controller has available noisy 
absolute position and attitude measurements that are used as input of 
the aircraft Kalman filter. The measurement (output) matrix is then; 


1 0 0 0 0 0 
0 1 0 0 0 0 
0 0 1 0 0 0 


(7.1) 


From data available in [26] and [27] , the intensity matrix of the measure- 
ments noise is selected to be: 


0 - 


0.286 0 

0 0. 001 0 

0 0 0 


•1 

n f 


(7.2) 


J 


With this choice of measurement noise, and with the driving noise 
resulting from the wind turbulence model, the Kalman filter poles are 
located at; -0.5743 ± j 0.5593, -1.317 ± j 1.104 and -1.107 ± j 1.318. 


- 151 - 



If measurements of the angular jities, as well as 

angles and disx^lacements, are availaLj ; navigation system, a 

Kalman filter can be designed, assuming the measurement of all the 6 
states. The measurement matrix is then the identity matrix: 

- I (7.3) 

The measurement noise intensity matrix is selected, using the same 
references as before, and is equal to: 

0.286 0 0 0 0 0 

0 0.001 0 0 0 0 

0 0 0.001 0 0 0 

000 0.286 0 0 
0000 0.00024 0 

0 0 0 0 0 0.00024 

In this case, the closed-loop poles turn out to be at: -0.4981, 

-0.5662 ± j 0.5804, -0.91, -3.92, -4.819. 

The first design is referred in the next sections as the 3 measurements 
case, and the second as the 6 measurements case. 

Note that in the situation when the 6 measurements are available, 
we are very close to the full state feedback. The role of the Kalman 
filter is not to reconstruct the unavailable states, as it was the case for 
the ship Kalman filter for example. In this case, its role is to provide 
optimal estimates of the states from noisy measurements, using the known 




- 152 - 



dynamics of the aircraft: it is actually nothing else than a filter. 

In the closed-loop system, its role is to filter the noise and so, it 
decreases the bandwidth of the system. The consequence is, theoretically 
(i.e. in the assumption of white measurement noise), an increased perfor- 
mance compared to the full (noisy) state feedback, at the expenses however 
of robustness, as will oe seen in the next chapter. 


7.4 LQG Controller Performance 

The evaluation of the controller performance can be done as for the 
LQ controller in chapter 6, by computing the rms values of the tracking 
errors and of the controls. 

To determine the rms values of the states (ship, aircraft, and 
Kalman filter states), we write the overall system equations: 


X 


• • 

X 


' i 


= A' . 


+ 


A 

X 


A 

X 


hQ. 


where; 



-BG 

A-BG-HC 


(7.5) 


(7.6) 


The overall state covariance matrix is obtained by solving the Lyapunov 
equation: 


A'^^P + PA' + 


0 


(7.7) 
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with: 



0 H0h'^ 


and: 

T 

E ) 
EiSix^) 


E ) 
E (^'^) 


(7.8) 


(7.9) 


The rms values of the tracking errors and controls are the square roots 
of the diagonal elements of: 


E (w, „w, 
'-A-S-A- 


J = W, 


A-S 


E(xx 


T. ,,T 


"A. 


(7.10) 


and: 


rp .. T T T 

E(cc ) = T G E (xx ) G T (7.11) 

— - c ~ c 

respectively. 

Tables 7.1 and 7.2 show the results obtained for the rms tracking 
errors and controls. In the case when the aircraft navigation system is 
of quality such that the navigation errors are negligible, the aircraft 
total state Kalman filter is not necessary. The ship Kalman filter is 
still required as a state reconstructor. The degradation in the tracking 
errors (compared with the unrealizable full ship and aircraft state 
feedback case) is then very small, while the rms controls required are 
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SEA 


H=12ft W »0.4807rad/s 
m 


Rms 

and 

Tracking Errors 
Controls 

Sway 

(ft) 

Roll 

(deg) 

Yaw 

(deg) 

w 

6t 

1,2 

(^) 

6a 

3 

(deg) 

bQ 

(1) 

0.801 

6.410 

0.329 

9.72 

_j 

16.70 

9.37 

LQG 

(2) 

0. 913 

6.432 

0. 341 



9.36 

LQG 

(3) 

2.312 

7.249 

2.54 7 

11.39 

23.61 

10.42 

LQG 

(4) 

1.401 

6.567 

1.499 

10.10 

18,13 

9.61 


(1),(2),(3), and (4) see below 

Table 7.1 ; LQG Controller Performance (decaying sea) 


SEA H=10ft 0) '=0.72rad/s 

m ' 


Rms 

and 

Tracking Errors 
1 Controls 


1 

Roll 

(deg) 

Yaw 

(deg) 

‘^'^1,2 

(ueg) 

6t 

° 1,2 
(%) 

6a ^ 
(deg) 

LQ 

(1) 

0.691 

2.480 

0.317 

4.72 

8.78 

3.92 

LQG 

(2) 

0.812 



4.70 
i 

8,76 

J 

3.89 

LQG 

(3) 

2.274 

4.193 

2.545 



6.01 

LQG 

(4) 

1.338 

2.856 

1.496 

5.45 

11.19 

4.47 
1 


(1) : LQ Controller, with wind disturbances 

(2) : LQG Controller, with only the ship Kalman filter 

(3) : LQG Controller, 3 measurements case 

(4) : LQG Controller, 6 measurements case 


Table 7. 2; LQG Controller Performance (fully developed sea) 
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sometimes smaller, due to the decrease in bandwidth caused by the ship 
Kalman filter. 

When the aircraft navigation system is of lower quality, and an 
aircraft total state Kalman filter is used, along with the ship Kalman 
filter, the degradation in porfomiance is more serious, but still 
reasonable in the b measurements case. 

7 . 5 Suninary 

In this chapter, we considered the optimal estimation of the aircraft 
states from noisy measurements. Two cases were considered: the first 
corresponded to measurements of position and attitude only, the second to 
measurements of position, attitude, and their derivatives (full state 
measurement). The second case is considered the nominal one. 

The degradation in performance due to the estimation of the ship 
states and of the aircraft states from noisy measurements was indicated. 

It was shown to be very small if only the ship states had to be estimated, 
and reasonable in the case of the noisy measurement of the full aircraft 
state vector. 
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CHAPTER 8 


ROBUSTNESS ANALYSIS 


8.1 Introduction 

The robustness of a feedback system to modelling errors (parameter 
uncertainty, luunodelled dynamics, neglected couplings ...) is certainly 
as important as the stability of the nominal closed-loop system. It is 
fundamental for aerospace applications in which, very often, the nominal 
open-loop system is unstable. While the open-loop system is an optimally 
robust system when it is stable, the existence of an optimally robust 
closed-loop system can be conceived in the case when the open-loop 
system is unstable. 

Although the issues of robustness are well understood, and easily 
expressed for single-input single-output systems, they are much more 
complex for multiple-input multiple -output systems. Much research has 
been done recently in this area, showing mainly the usefulness of 
singular values to quantify robustness properties of multiple-input 
multiple-output feedback ;.ystems. In this thesis, we will mainly refer 
to the results found in [23] , [24] , [28]. 

8.2 Robustness Measures 


8.2.1 Sin3^e2input_S^ng_le-Ou^ut_ 

The stability of SISO feedback systems can be determined by the use 

of the Nyquist diagram. The Nyquist contour D is defined as two segments 

JK 
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of the positive and negative imaginary axis, connected by a half circle 
whose radius R is, to the limit, taken to be infinite, so that the whole 
right half plane is then included in the contour. The Nyquist contour 
transformed by the complex function g (s) (the open-loop transfer function) 
is then plotted in the complex plane. The closed-loop system is guaranteed 
to be stable if the number of counterclockwise encirclements of the (-1,0) 
point is equal to the number of open-loop unstable poles. 

Robustness is easily measured by the distance from the transformed 
contour to the (-1,0) point. It is usually expressed in gain and phase 
margins. They indicate what minimal multiplicative perturbation (constant 
gain or phase shift) would make the Nyquist diagram pass through the 
(-1,0) point. This would be at the limit to change the number of encircle- 
ments, and, consequently, the stability of the closed-loop system. 

8.2.2 Muiti^^e2^Ir^]£t_Mul^i£^l£-Ou^ut_(MI5ra Ga&e_ 

The stability of MIMO feedback systems can be determined in a very 
similar manner by a Nyquist diagram. The transfer function to be plotted 
on the Nyquist diagram is now (-1 + det(I+G(s))) where G(s) is the open- 
loop transfer matrix. The same criterion can then be used. 

We note immediately however that det (I + kG(s)) is different from 
1 tkdet (G (s) ) . This indicates already that no stability margins can be 
found from the MIMO Nyquist diagram, as in the SISO case. 

Very often, separate loop-by-loop stability margins, similar to the 
ones used in SISO systems, are considered indicative of the overall system 
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robustness. They are however clearly insufficient. Perturbations are not 
likely to occur separately, and simultaneous perturbations in different 
channels may be much more dangerous for the stability than separate 
perturbations. The importance of possible cross-channel disturbances has 
also to be evaluated, and is not reflected by a single loop-by-loop 
robustness analysis. Consequently, going back to the SISO case, with 
a single loop-by-loop analysis is not a satisfactory way to study the 
robustness of MIMO feedback systems. 

Unfortunately, the distance from the (-1,0) point to the Nyquist 
contour transformed by (-1+det (I+G (s) ) ) - which is the same as the 
distance from det(I+G(s)) to 0 - is not a reliable measure of robustness. 
In fact, numerical analysts have long recognizied that the specific value 
of the determinant of a matrix is a poor measure of the nearness of this 
matrix to singularity (or rank deficiency). In other words, det(I+G(s)) 
can be large, although a small additive perturbation E can make 
det(I+G(s)+E) equal to Zero. 

8.2.3 Singular _Value£ 

These shortcomings can be avoided by the use of singular values. The 

singular values of a complex matrix A are defined as the square roots of 

H H 

the eigenvalues of the matrix A A (A is the complex conjugate transpose 
of A) , i.e. : 

0^(A) = xJ(a“a) (8-1) 
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Th« maximum singular valui^ of a matrix A is indicative of its "size"/ 


it is equal to the spectral norm of the m'^tX'ix A: 

II Ax |L 




max 


( 8 . 2 ) 


3^0 


where : 

II X II 2 o (8.3) 

The minimum singular value of a matrix A is indicative of its nearness 
to singularity. If 0 , -QL, there exists a perturlation E such that o „„(E)*a 

mi.n uieix 

and det(A+E)=0, Moreover, no matrix E such that 0 (E)<a will maHe A+E 

BioX 

singular. The minimum singular value is then an indication of the minimal 
"size” of the additive perturbation required to bring the niatrix A to 
singularity. 

An advantage of the use of singular values is that it also allows to 

compute the minimal perturbation E through the singular values decomposition 

H H 

If U and V are matrices containing the unit eigenvectors of AA and A A 

(known as left and right singular vectors of the matrix A) , and if £ is a 
diagonal matrix containing the singular values of A, then A can be 
decomposed as: 

A = U Z v“ (8.4) 
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And tht matrix B can be decomposed as: 


**8 H 

. V" (8.5) 

0 *•0 I 

nj 

where P_ is any n-1 x n-1 matrix such that || P_ II . < a and where it is 
s ' » " 2 n 

assumed that the singular values have been ordered so that a > a . (A) . 

n min 

Without loss of generality, one can select P-»0, so that: 

s 



8.2.4 Spbustnes s_o£ 

The robustness of MIMO systems is defined by the size of the minimal 
perturbation E that brings the closed-loop syst«n at the limit of instabi- 
lity. Then, the return difference matrix ItG(e) (where G(s) is the 
perturbed open-loop transfer matrix) is singular for some s on the 
Ny«:juist contour. 

Different structures of perturbation can be imagined, for example, an 
additive perturbation: 


8 - 6 -f E 


(8.7) 


or a multiplicative peirturbation: 
8 - G (1-tE) « G.L 


( 8 . 8 ) 
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In tho first case, the closed-loop system stability is guaranteed if: 




(6.9) 


In the second case, if: 




> 0 (E) 

max 


( 8 . 10 ) 


Other structure of errors have been analyzed in (23 ] and showed the use 
of ( (ItG) (1-G) for some type of j)erturbations. 

As in th( SISO case, the transfer matrices go to zero as s-+<» , so 
that the only important of the Nyquist contour is the imaginary axis, 
and s can be replaced by ju) in the previous formulas. 


8.2.5 Gu^ran^O£d_Ga^in and jPhas^ 

Guaranteed gain and phase margins can also be obtained from the 
minimum singular values. Although only diagonal perturbations are 
considered here, they are not limited to the case when only one channel 
is perturbed at a time. At the contrary, simultaneous gain or phase 
changes in all channels together are considered. The results presented 
here can be found in [23] . The gain and phase margins obtained are 
guaranteed (or conservative) margins. In other words, the gains (or phases) 
in all channels can be changed simultaneously within the limits indicated 
by the gain (phase) margins without causing instability, but it is not 
necessarily possible to find some perturbation of this amplitude that 
will cause instability (at the contrary of the previous case where a 
general type of perturbation was considered) . 
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For the different type of errore considered in [23 1, the guaranteed 


gain and phase margins are: 







(8.11) 

PM o <-2sin^ 1 , 

-1 

ot. 

(8.12) 

2sin 

2' 

-1 




Oj^in(i+G ) > a GM 3 (1-a , 1+a) 



(8.13) 

PM 3 (-2sin^ j , 

2sin^ 

-) 

(8.14) 

Omin((I-«5) (X-G)-^) > o -*■ GM 3 (~ 

14d 
’ l-o' 

(8.15) 

PM 3 (-2tan^a , 

2tan^a) 

(8.16) 

8.3 VTOL Control System Robustness 




8.3.1 Introduction 





The robustness results sunstiarized above are applied to the LQ and LQG 
control system designs of chapter 6 and chapter 7. Note that the computa- 
tion of singular values of transfer matrices can be done quite rapidly, 
by the use of efficient algorithms [29], [30] , [ 31L 
For the LQ design, we have: 

X » Ax + Bu (8.17) 

u = -Gjc (8.18) 
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The opan-lcx>p transfer matrix isc 


G(jO») - G B (8*19) 

We can also consider the transfer matrix: 

G' (jt^) » R** G(jW) r"** (8.20) 

which is equivalent to consider a normalized transfer matrix (normalized 
by the weights attributed to each control in the F matrix). 

The Kalman inequality is: 


(I+G' (jW)“ (I-KS'(jW))^ I (8.21) 

This guarantees that: 

Omin(i+G' (jW)) ^ 1 (8.22) 

and then: 

GM D (1/2, “») (8.23) 

PM D (-60deg,60deg) (8.24) 


As: 

+1 > C , (I-K3'"^) (8.25) 

rain — min 

This also implies that: 

(I+G'“^(jU))) > 1/2 (8.26) 
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For the LQG design, we have: 


• 

X - 

Ax + Bu 

(8.27) 



U B 

-G^ 

(8.28) 



' • 

X " 

(A-BG)x + H(^-Cx) 

(8.28) 

y a 

Cx 

(8.30) 


so that: 


G(jW) - (0 G). 



r 


-1 

• ■ 


A 0 



B 

jO)I - 

HC A-BG-HC 


• 

0 





L 


(8.31) 


Again, to have some normalization of the transfer matrix, the matrix: 


G’ ( ju») 


= R** G(jW) r"** 


(8.32) 


can be used. In the LQG case, there is no inequality as (8.22), (8.26), 
and, consequently, no guarantee on the robustness of the closed-loop 
system. 

In the next subsections, the minimum singular values of the loop 
transfer matrices for the and LQG designs are indicated. To provide 
some normalization, the loop transfer matrices considered are those given 
by (8.20) , (8. 32) (note that the parameter p has no influence in this 
normalization) . 
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8.3.2 |42_MiniJnuin_Sinc|ul^a£ Val^ues_ 


The minimum singular values of the LQ desisin are shown in Figures 
8.1, 8.2, and 8.3. It is easily checked that Oj^^^(I+G ^)^l/2 and 


a . 
min 


(I+G) > 1. 


It is interesting to note that the minimum over U> of the minimum 
singular value (for any perturbation criterion) is smaller for p= 30 
than for p= 0.3. In other words, the closed- loop system with higher 
loop gains is more robust than the others. This is probably due to the 
instability of. the open-loop system. 


8.3.3 yjG Minimum £ingular_ V£lue£ 

The minimum singular values of the LQG design with 3 measurements 
(aircraft position and attitude) are illustrated in Figures 8.4, 8.5, 
and 8.6. For the 6 measurements case (motions and derivatives), the 
values are shown in Figures 8.7, 8.8, and 8.9. 

The first conclusion is that the robustness is seriously degraded 
when only position measurements are available. The estimation of positions 
and velocities from noisy measurements of positions only does not only 
result in degraded performance, as indicated previously (tables 7.1 and 
7.2), but also in reduced stability margins. 

The robustness is improved when the 6 measurements are used in the 
Kalman filter. However, the minimum singular values drop significantly 
for very low frequency in the nominal case p=0.3. 
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Figure 8.1 ; LQ Minimum Singular Values 
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Figure 8.2 ; LQ Minimum Singular Values 
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Figure 8.3 ; LQ Minimum Singular Values 
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LQG Minimum Singular Values (3 measurements case) 
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8.7: LQG Minimum Singular Values (6 measurements case) 
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8-8; LQG Minimum Singular Values (6 measurements case) 
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The singular value decomposition can be used to obtain the minimal 
destabilizing perturbation, and get some insight about what physical 
effects ate implied by this perturbation. The minimal destabilizing 
perturbation, computed for p^O.! and u)bO, depends on the error criterion 
selected, but it turns out that the different results are very similar. 
The minimal multiplicative perturbation, indicated by ^>, is; 



0.089 

0.063 

-0. 065 

tl « 

0.0785 

0.955 

0. 046 


0.133 

-0. 076 

1.079_ 


This can be interpreted in many different ways. For example, we can 

assume tliat this perturbation i.s due to modelling errors in the B matrix 

(force generation part, or control effectiveness). This would mean a 

change from B to i =B L. The nominal B was given in Table 5.1; 
c c c c 


B * 
c 



0 

0 

0 

0 

0 

0 

0 

0 

0 

21.231 

0 

-10.989 

2.976 

4.595 

1.265 

2.864 

0.194 

-2.685 


(8.34) 
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The matrix B is then: 
c 


B - 

c 


0 

0 

0 

•20.334 

3.438 

2.217 


0 

0 

0 

•0.290 

4.443 

0.410 


0 

0 

0 

-10.468 

1.536 

-3,071 


Which, for the last three rows, is a change of: 




-4% 

CO 

-5% 

% , . - B , 





> c,i] c,i; 

c,i] 

^ 

16% 

-3% 

21% 



-23% 

111% 

14% 


(8.35) 


(8.36) 


We see that the "destabilization" requires important changes in terms that 
were quite small originally, and could be called parasitic. For example, the 
second term of the last row is increased by about 100%. We remember that 
this term expresses the yaw acceleration produced by a roll command, and 
that this is due, originally, to the angular difference between the body 
axes and the principal axes. 

At this point, the designer will decide, on the basis of the knowledge 
of the accuracy of the model, whether such perturbation (or error in the 
model) is physically possible or not. This is a delicate task tf:at requires 
a good knowledge of the systenv and of the precision to which the model is 
known. In our case, no data could be found in [ 1] , and in the original 
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model (19l about the precision of the model, so that this question could 
not be answered with certitude. However, small singular values over some 
frequency range are indicative of decreased stability margins of the closed- 
loop system, and the best solution is, if possible, to clkange this charac- 
teristic: it is the subject of the next section. 

8 , 4 .Robustness Re c overy 

The robustness recovery for bQG designs refers to he recovery of the 
LQ transfer matrix - and consequently stability margins - from the LQG 
design. The robustness recovery technique was proposed in f28), following 
a procedure very similar to the one proposed in [32]. 

The idea is simply to modify the Kalman filter in order to approach 
the liQ design. The filter is redesigned with a new value of the driving 
noise spectral intensity matrix; 

S ® E + q (8.37) 

o 

This corresponds to assume that some white noise is present at the input, 
and, consequently, tells to the mathematical expression of the problem 
that there are uncertainties at this point. It can be showni 28] that, as 
q-w, the Lg loop transfer matrix is asymptotically approached by the LQG 
loop transfer matrix. 

The robustification procedure is applied to the LQG design described 
in chapter 8 (6 measurements case), and is illustrated in Figures 8.10, 

8.11, and 8.12. For q«0.01, the robustness properties are seriously 
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Figure 8.10 ; LQG «laimum Singular Values (6 meas. , with robustif ication) 
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improved, and for q°'l the LQ singular values are almost recovered. 


The minimal values of the minimum singular values over the whole range 
of frequencies, and the gain ^md phase margins obtained from equations (8.11) 
to (6.16) are shown in Table 8,1. As these margins are guaranteed, or 
conservative, the overall guaranteed margins are the union of the margins 
corresponding to the different criterions, and are underligned in Table B.l. 
For q»0.01, the guaranteed gain margin is about from 0.4 to 3, and the 
guaranteed phase margin about ±40deg. These appear sufficient for our 
application. 

The robustness recovery is however made at the expense of the perform- 
ance. In fact, the bandwidth of the Kalman filter increases quickly as q 
increases. Table 8.2 indicates the location of the Kalman filter poles for 
different values of q. 


q 

Aircraft 

Kalman 

filter poles 



0 

-4.82, 

-3.92, 

-0.91, 

-0.57 i 

j 0.58, 

-0.5 

0.01 

-10.81, 

-5.37, 

-2.23, 

-1.12, 

-0.60, 

-0.49 

1 

-97. , 

-43. , 

-7.49, 

-1. , 

-0.5 , 

-0.49 

100 

-989. , - 

429. , 

-69. , 

-1. , 

-0.49, 

-0.49 


Table 8.2 ; Aircraft Kalman Filter Poles (with robustif ication) 
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The size of the loop transfer matrix G is represented by its maximal 
singular value, so that the bandwidth of the contro.' system can be defined 
by the crossover frequency of the transfer function ( jW) ) . The 

evolution of this transfer function with the parameter q is shown in 
Figure 8.13. The increase in bandwidth resulting from the robustif ication 
procedure can be directly observed from this figure c 

The increase in bandwidth may result in improved or degraded tracking 
performance, but almost certainly in increased control authority, as more 
noise is passed through the filter. Table 8.3 illustrates the changes in 
the rms tracking errors and control authority for different values of q. 


SEA H=10ft 0) =0.72 

m 

rad/s 





Rms Errors and 

Sway 

Roll 


6 o( 

“l ,2 

■ ^ 1,2 

603 

Controls 

(ft) 

(deg) 

(deg) 

(deg) 

(%) 

(deg) 

LQ w.A/C noise 

0.691 

2.480 

0.317 

4.719 

8.78 

3.917 

L 2 G q =0 

1.338 

2.856 

1.491 

5.448 

11.19 

4.467 

q= 0.01 

1.350 

2.869 

1.482 

5.456 

11.83 

4.653 

q=l 

1.239 

2.851 

1.445 

6.394 

18.70 

5.825 


Table 8.3; LQG Controller Performance (6 meas., with robustification) 
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For q“0.01, the decrease in performance is very small, while for q=l 
it becomes significant. The design corresponding to q=0.01 appears to be 
a very satisfactory one, both for tb; performance and the robustness. 

Note however that in this performance evaluation, white measurement noise 
has been assumed. In reality, the noise will be bandlimited, and the 
degradation will probably be less than indicated. In fact, as all the 
states are available (but noisy) , an LQ controller, without Kalman filter 
(using the known dynamics of the aircraft), is not an unrealistic possi- 
bility: its pf rformance would be less t)ran the performance of an LQG 
design, but the robustness would be increased, and the computational 
load would be dramatically reduced (this would be the case if a good 
navigation system is available aboard the aircraft, cf. chapter 7). 

8.5 Nyquist Diagrams 

Some Nyquist diagrams will illustrate the facts indicated at the 
beginning of this chapter. Figure 8.14 shows the function -1+det (I+G ( jw) ) , 
0)^0. To plot it, it is convenient to compress the distances radially by 
a logarithmic transformation which is chosen to be: 

r = log (1-f/ x^-fy^' ) 

x' = X. r 

y' = y.r 


(8.38) 

(8.39) 

(8.40) 
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For MXMO systems « such compression of the distances is almost 
unavoidable, as the transfer function det(X4-G(s)) is of an order equal 
to the order of the system (equal to 6 here) . 

Figure 8.15 shows the complete contour, and the number of encircle- 
ments for the LQ design. There are 3 counterclockwise encirclements, and 
1 clockwise encirclement. Remembering that there are two unstable open- 
loop poles, this confirms the closed-loop stability of the. LQ controller. 

As in the SISO case, some attention has to be given to the two poles at the 
origin, and to the encirclement at the infinity. 

Figure 8.16 shows (in normal cartesian coordinates) the Nyquist 
diagrams for the LQ and LQG design (6 meas. ,q=«0) near the critical point. 

It is checked that the Nyquist diagram for the LQ design never enters the 
unit circle centered at -1. This is a consequence of the Kalman inequality 
(8.21) [33] . 

The minimum singular values shown in the previous subsection indicate 
the size of the minimal perturbation that destabilizes the closed-loop 
system. It is also possible to compute the value of the multiplicative 
(or additive) perturbation L, such that the perturbed transfer matrix G.L 
is at the limit of instability. In terms of Nyquist diagrams, the minimal 
perturbation L is the perturbation that , at the specific frequency is where 
it is computed, brings the point of the Nyquist diagram -1+det (ItG ( jw) ) 
to the -1 point (then det (I+G ( jiS) )=0 ). At any frequency, the minimum 
singular value indicates the size of the minimal perturbation that will 
make the Nyquist diagram pass through the -1 point, and bring the system 
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at limit of instability, over the whole range of frequencies* the 
minimal value of the minimum singular value gives* from all the minimal 
perturbations found at each frequency* the constant perturbation chat has 
the minimal size. The meaning of the previous discussion* is that the 
frequency at which this overall minimum occurs may correspond to a fxjint 
located very far from the -1 point on the Nyquist diagram. This is 
illustrated in Figures 6.17 and 8.18. 

Fig. 8.17 show the Nyquist diagram of the loop transfer matrix G 
corresponding to the robustified (q=0.01) LQG (curve LQG) together with 
the diagram of the loop transfer matrix G.L (curve perturbed LQG) . The 
perturbation L is the minimal multiplicative perturbation that corresponds 
to Pig.6. ' r.. It is computed at the frequency where the minimum singular 
value in minimal* i.e, 1.15rad/s. The distortion of the Nyquist diagram 
is quite peculiar. The point that is brought to the -1 point is not at 
all the closest point of the Nyquist diagram* although it is the closest 
point on the basis of minimal multiplicative perturbation. This is a 
practical proof of the fact that the distance to the -1 jxaint in the 
Nyquist diagram is not representative of the closeness to inKtability, 
with respect to general type of perturbations. 

The situation is even more surprising for the non-robustified L^^G 
design (q=0) . In this case* the minimum singular value goes to its minimum 
at o)=0 (cf. Pig. 8. 7). At this frequency* due to the presence of the two 
poles at the origin, the determinant of I+G(ju)) is infinite. This means 
thiTit the closest point (closest in the sense of minimal norm additive or 


- 191 - 



or POOH 


wuai. jry 



Figure 8. 17 : LQG and Perturbed LQG Nyquist Diagrams, Robustified 
Case with q=0.01 (diagram radially compressed) 
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muK'iplicatlve perturbation of G(j(i})) from the Nyqulcit <Ha<Jr^un to the -I 
point i> the s>oint at inflnityl This .illustrated in Fig. 8. 18. As this 
vase i'3 numerically perfectly ill-conditioned, the minimal multiplicative 
perturbation is computed at u)*0.1 rad/s, where the minimum singular value 
is very close to the minimum value reached at a)>»0rad/8. 

8.6 Summary 

In this chapter, we addressed the important issues of the control 
system robustness to modelling errors and parameter uncertainty. Due to the 
decoupling between the ship and the aircraft, this problem only involved 
the aircraft model and the feedback from the aircraft states. 

The importance of the singular values analysis in the definition of 
multivariable robustness measures was first introduced, and rnme recent 
results concerning the stability margins of multivariable systems were 
summarized. 

The singular values analysis was used for the different transfer 
matrices of interest, and for different designs considered previously. 

For the ncaninal design, it was shown that reduced stability margins were 
obtained at very low frequency, and some physical interpretations were 
indicated for the minimal multiplicative destabilizing error given by the 
singular value decomposition. A robustification procedure was used, and 
appeared successful in retrieving the favorable robustness properties 
of the LQ design, with a limited degradation in performance. 
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Finally, some Nyquist diagrams illustrated the theoretical results 
mentioned previously, especially the impossibility of defining stability 
margins from the Nyquist diagram in the MIMO case, in contradistinction 
to the SISO case. 
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CHAPTER 9 
CONCLUSIONS 


9. 1 General Conclusions 

Two major aspects of the problem of landing VTOL aircrafts on small 
ships were addressed in this thesis. The first is the estimation and 
prediction of ship motions, mostly for applications in piloted landings, 
and the second is the design and evaluation of an automatic controller, 
for use in fully or partially automatic landing schemes. 

An accurate ship model was first der 'red, using hydrodynamic data 
Cor the DD963 destroyer. It was shown that the ship motions hcive power 
spectra that are concentrated in a narrow frequency band, especially the 
roll motion. Consequently, good prediction time could be obtained, and 
the roll motion could be predicted for as far as 5 to 10 seconds ahead. 
This was a significant result, as the roll motioif turns out to be the 
most important lateral motion in the VTOL landing problem. In particular, 
the lateral translation of the landing pad is due to a large extent to 
the roll motion, as the location of the landing pad is significantly 
higher than the ship center of rotation. 

The sea modal frequency is a significa»;t parameter of the ship model, 
and its large influence on the estimation and prediction errors justifies 
a precise estimation of this parameter in a real-time application. 

The aircraft model was briefly described, and the important couplings 
amongst the motions were mentioned. Two of them appeared significant, and 
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contradictory: the first is the lateral acceleration due to a roll angle 
of the aircraft, and the second is ^,ho roll moment due to a lateral 
deflection of the thrusts. 

The design of an automatic controller, defined as a tracker of the 
ship motions, made clear some other of the constraints related to the 
tracking of the ship motions. The relations between the lateral transla- 
tion of the ship landing pad induced by the roll motion at the center of 
rotation, and the roll motion itself, appeared to be contradictory with 
the two important couplings of these motions in the aircraft itself. 

These constraints are specific to the landing problem, and are independent 
of the control system design metliodology . Although the roll motion may 
often be neglected, its possibly large amplitude in high sea states makes 
imperative to take these limitations into account in an engineering design. 

The design of tlM3 automatic controller also illustrated some results 
and recent advances in the design of an optimal controller, using the LQG 
methodology. The relations between the choice of the state and control 
weighting matrices in the quadratic cost, and the resulting root-loci ware 
explained, and justified the careful selection of these matrices. 

The usefulness of the singular values analysis in the definition of 
realistic stability margins was indicated, and some Nyquist diagrams 
illustrated related theoretical results, and the important differences 
between the MIMO and the SISO case. 

The robust if ication procedure appeared to be very successful in the 
design of robust LQG control systems, and led to increased stability 
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margins, with moderate degradation in the i>errormanco. 


Research 

Among the .possible topics for future research, we note: 

- the evaluation of the performance and of the stability margins of 

an automatic controller including the longitudinal and the lateral motions, 
and the study of the influence of the cross-coupling terms (gyroscopic 
terms ) 

- the inclusion of the aerodynamic effects in the aircraft model, 
and of the ship airwake turbulence and ground effects; the aerodynamic 
effects being nonlinear, some nonlinear theories may then be useful 

- although this research concentrated on the landing problem, some 
interesting results may be obtained from the design of a control system 
for the transitiot: flight, and from the use of modern control theories 


for this time-varying problem. 
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APPENDIX A 


A. 1 Ship Model Values 

A. 1.1 Sea 

J ■=• 0.707 


U X 

a a= — 0» CX)8® 

g m 


SHIP MODEL COMPLEMENTS 


B(a) 


1.9339 
14- 2a 


S • 0.3125 7 “ B (a) 
o 0) 

m 


u 


Wjj|(l+a) 


b 0.8409 


A. 1.2 Porce Dyneun^ics_ 
UbO, (|)=590 degrees: 


F “ 310 
s,o 


F <= 2120 
r,o 


F = 11300 
y,o 


= 0.72 

U) = 

0.6 

s,o 

s,o 


„ = 0.7 

(i) = 

0.76 

r,0 

r,o 


= 0.35 

0) =5 

0.96 

y»o 

y,o 



Any U, (j): 


F » P Sind), J » J _ sind), 0)„ = 4 0)^ „ ^ cosb ) sincj) 

s 8,0 ^ s s,o ^ s s,o s,o g ^ ^ 

F^ = F^^^ sin<J), Jj. = sin(j), 0)^ “ ( cos(^ ) sine)) 

” *'y,o *^y " ‘^y,o “y “ ^ ‘“y,o “^y,® | 
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A. 3. 1 S^hy3_pjfnam^c8 


U«o, 4>“90 deyrees; 


215 

988 

-230 

908 

104000 

0 

-230 

0 

3.76 10 

223 

-759 

14600 

-759 

22900 

182000 

14600 

182000 

O 

CO 

10.6 

-55.4 

423 

-55.4 

887XW3P 

6270 

423 

6270 

144000 

C 

0 

0 

0 

28800 

0 

U 

0 

0 


RDF=3 (the roll dampiny factor takes the nonlinear dampiny into 
account) 


Any U , 4’ ! 


A. 


‘h,o 

P 


0 

0 


0 

0 


h,o 


1,1 


h,o 


-B, -B. U.A. 


1,2 

1,1 
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Bv- “ B. + U . 
h h,o 


0 

0 


0 

0 


■■A, 


h«o 


1,1 


h,o 


1,2 

A, A. U B, 

p 


where : 


U COBd) 

0) « 0.425 + 

p 178.27 


A. 2 Ship Dynamics; State-Space Model 
We have the equation of motion: 


* »lA. * =h2h ' 


where: 


V . 

--h 


sway 

roil 

yaw 


F = 


F 


sway force 
roll moment 
yaw moment 


We call: 


R = (M->-A. = {r. ,} 

h 13 
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P • (M+A^^) «* tp^j) 

Q « (M-^A^ » R.C, 

h h h 


Roil ia the only motion havinq a sprinq constant so that: 




0 

0 

0 


rr 

0 


0 

0 

0 


and: 


Q 


r, „c 0 

12 rt 

r„„C 0 

22 rr 

r ,, 0 

j 2 ri. 


Uslny Laplace transforms, we write: 


or: 


s"x = 

-Psx - Qx + ; 

s (x^ ) 

= (SXj^) 

S(X2) 

= (sx^) 

S(X3> 

= (SX3) 

S {sx^) 



«■ 


12 


12 rr '^^2' 


P„(s«3> + 


(A.l) 
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(A. 2) 




8(SXj) 


- - *^22 "*22‘^rr^*2^ " *^23 '■*3^ 


^ *^2/1 '*' *^22*^2 * *’23^3 
- ‘?31<®’'l> " P32<®^2> ■■"32''rr<’*; 


■*^ *^31^1 ^ ‘^32®’2 ■*■ *^33^3 


^33 ^***3^ 


(A. 3) 


which is equivalent to a state-space representation having 6 states. 
There is an implicit pole-zero cancellation at thai origin, between the 
the poles ot the sway and the yaw motions and the zeroes of the corres- 
ponding force and moment. 

To avoid this, we rewrite eq. (A.2): 


IK 


Pr 


rr s 


1 *^21 *"22 *'23 

■ I M U m i ■ <3V ■■ jm — I. — ! I. V ■■ mrnr^nmntt V M» «» ■—■■■'— V 

^22 *■ *'22 ^ 22 ^ 22 


, *^21 ^1 , 2 . *^23 ‘''3 

-22 « 


(A. 4) 


and use (A. 4) with (A.l) to obtain: 


”<* 1 ' " ti2(«j>+ ‘ijlXj) ♦ 

" “U'T' * “15 


(A. 5) 


where ; 


11 


'12 


“13 


12 

-22 

-12 

-22 

-12 

r *^22 
22 


- P 


11 


- P 


12 
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tl iM 

14 

'12 

■■22 

*^13 


U - ■* 


r , 


11 

hi - “r^^ 

' 


U n 

15 

r - ^ 

“ -^22 

*"23 


Similarly, 

(A. 4) and 

(A. 3) 

give; 

s(X3> 


^ <^42 

<«2> 


'44' 3' 


^41^ s ^ *'■ “45 ^'s ^ 


where : 


*^41 “ ^21 " ^3X 


t . _ = 


•32 


42 r 


22 


32 

"43 " ^22 - P32 

. ^32 

"^44 “ P23 " ^^33 

^32 

n 8B r — ■ "" ’ r 

41 *31 r,.^ ‘21 

22 

*^32 

^45 " ^33 “ r^2 *^23 


And, of course; 


(A. 6) 


S^X^) (sx^) 


(A. 7) 



rinally, (A. 4), (A. 5), (A. 6), and (A. 2) lead tot 




* “22 ♦ “24<^2> * “25<T> * “26‘''3> 


(A. 8) 


where: 


"21 “ - PaHi - ^23^41 


^22 “ ” ** 21^2 ‘ ^ 23^^42 " ^22 


^23 “ - ^2l'^13 " ^23^43 ' *^22^rr 


♦^24 ‘ " 


21 1-1 *"23 44 


“21 “ “ ^21 '^11 " ^23^41 


^22 * ^21 


U Bt IT 

24 22 


^^25 " " ^21^15 ■ ^2S'^45 


u r 
26 23 


Equations (A. 4), (A. 5), (A. 6), (h.1) , and (A. 8) constitute the new state- 
space representation f with only 4 states: 
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’‘l 


’'l 


P/s, 





p. 

sx. 




1 


“ T . 


+ U . 

*"2 



X„ 



2 


2 


p./s 





P. 

x^ 




3 

3 


3 




and: 

T = ft . , } 

U = fu, , ^ 
are given above. 

The pole-zero cancellations are eliminated by the introduction of 
the integrals of the sway force and yaw moment, which, instead of being 
computed (implicitely) in the ship dynamics model, can directly be obtained 
"drift-free" from the transfer functions between the sea height and the 
generated forces. 
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APPENDIX B 


DISCRETIZATION PJ«)BLEMS 


A sin\ple derivation is presented here that explains some dit'ficiiitiea 
that can be encountered in discretizing systems containing lightly damped 
second-order modes. 

Assuming the system: 

^ * A X 

having a second-order mcxie X(A) w -o i juii. The exact discretization of the 
system is: 

x(t+At) = A^ £(t) 

whore ; 

AAt 

A , «• O 
d 

It would lead to the equivalent mode in the z-plane: 

(-Oijui)At 

X(A^^) » e ' 

If we use an approxirtote representation: 
x(t+At) » (ItAAt)x(t) 
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the nodes will, instead, be at: 

Ml+AAt) - (1-oAt) ± jwAt 

Let us consider consider a system A', having inodes at -'a ± jb in the 
continuous time domain which, discretized, would have the same modes as 
1+AAt. 

The following equations can then be obtained: 

^-aAt _ ^jbit ^ 


and: 


In / (i-aAt)2 + (a)At)2' 


1 , -1, wAt , 

At 


It is easily checked that for At vanishingly small, a approaches a, and 
b approaches 0) so that the discretization is valid. 

The question is to determine how small At actually has to be. 

For b to approach O), we need; 

wAt << 1 and oAt « 1 


or: 


At « l/o) and At « l/o 


which is the expected condition that At must be small compared to the time 


cunutants ot. the system. 


The condition for a is a little different. We need: 
/ l-2oAt+(oAt)‘* + (u)At)^* / l-2oAt' a i-oAt 


so that the conditions are: 

2oAt « 1 , (oAt)^ « 2oAt , and (o)At)^ << 2oAt 


or: 


At « l/2o 
At « 2/0). (0/0)) 

The first condition is similar to the previous one, while the second is 
similar only if o/o) is close to one. If a/o) is much le^ss than one (lightly 
damped mode) , the required time step is much smaller than the one indicated 
by the time constants of the system. If this is not taken into account, the 
term {u)At)^ will influence the 2oAt term, which is equivalent to reducing 
the value of a. This implies a reduced value of the real part, while the 
imaginary part remains constant. In conclusion, the discretization of very 
lightly damped second-order modes with approximate discretization may lead 
to even less damped equivalent discrete modes, even if small time steps - 
compared to the time constants of the system - are used. 


APPENDIX C 


TURBULENCE MODEL 


A wind dialurbanc*' nKjdcl is avaiiablo in ll.9j and is sununarized here. 
The moan wind component: has been noqlccted (only the turbulence is 
cot>siderud) , and the air turbiilence considered does not include the ship 
uirwakc! and the ground eiCect.s. 

The turbulence is modelled as the outputs of Dryden filters driven 
by white noises of inteusity 1. The outputs are additional turbulent 
volocit.ios of the aircraft v_., (sway), p_., (roil), and r,,., (yaw). The 
exi ressions give.n in (19]are: 


p„ -y O' 
^2 V 


fU ' 1 + /? 77- S 

V V 

V L ~ 

( 1 + s ) 


BN 


V 0 


W 


li V 
w 


, W. 1/& 


1 + 


irv’ 


V 


BN 


BN 


•y - 


■' * w = 


^BN 


whore p,, and p 


are white noises of intensity 1 . 



In the condition considered in this thesis: 
h=“40ft (altitude of the stationkeeping point) 

V==15.5Ct/s (-^lOkts) 

(-20kts) 

And the parameters of the filters given in jigj are: 

L “600ft 

V 

0 ^“6 ft/s 

b“44.43ft 

0^=3.O58ft/s 

L “69. 1ft 
w 

TO bring this wind turbulence model to a standard state-space formulation, 
the outputs of the Dryden filters are first multiplied by s, and entered 
in the acceleration equation of the aircraft model, instead of the velocity 
equation. Next, the time constants of the shaping filters being large 
compared to the time constants of interest, they are neglected, so that 
the Dryden filters become simple gains; 



V “>• r 

BN 3b BN 
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And the driving noise spectral intensity for the aircraft model is tlien: 
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APPENDIX D 



ROW 

COL 1 

' 2 

3 

i 

-4.61340-01 

1.44580+01 

2,77020+00 

2 

-2.89140+00 

-5.57740+01 

-8, 364 10+00 

3 

1 ,04080+00 

3,97390+01 

6,44460+00 

4 

-3.83380+00 

-6,59960+01 

-9. 19900+00 

5 

3.77480+00 

6.59390+01 

8.431 10+00 

6 

-3.30590+00 

-5,02900+01 

-6,23640+00 

7 

8,11 170 + 01 

3.21810+03 

5.31510+02 

e 

- 1 , 39680+02 

- 1 . 74530+03 

-2,398 10+02 

9 

3,67230+02 

2,67460+04 

4,85040+03 

ID 

- 1 . 78170+03 

-2 , 86820+04 

-4,30070+03 

1 1 

- 1 , 55940+03 

1.78220+05 

4.087 10+04 

12 

-2. 18620+04 

-4. 16170+05 

-7. 16680+04 

13 

5.83240-01 

6,46560+00 

7.91960-01 

14 

- 1 ,04930-03 

1.57160-01 

3, 1 1490-02 

15 

1 , 29310-02 

4 , 77970-01 

4.73700-02 

16 

1 . 58390-03 

4.73700-02 

9 . 16890-03 
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Feedforward Ga^n_M£t£ix 

COU 1 23456 

ROW 

t !5.)087D-04 7.30470-04 l.6527D'04 -2.73530-04 5.11660-05 -9,90030-06 

2 1.74390-03 -1.90280-03 6.07390-04 -7,24010-04 2,07200-04 -2.77170-04 

3 3.45780-05 -1.29550-05 1,37760-05 -1.16100-05 4,91620-06 -6,17380-06 

7 a 9 10 11 12 

1.11080-04 4,77910-05 1,13460-06 -4,73370-06-7.57810-07 -2,49620-07 

3,00830-04 1,66790-04 4.38430-06 - 1 , 49240 -05 “ 1 . 88 190 -06 -8,00590-07 

1.04750-07 1,69060-06 6.99250-00 -1,63190-07 4.80160-08 5.34880-0? 

13 14 15 16 

4.54870-02 -1,04920*00 -6,39630-01 -5,40860+00 
9,28240-02 -3.75420100 -4,17850+00 -9.95670*00 
-1,09120-03 -2,62430-02 -1,41860-03 9,89060-01 


^ Aircraft Gain Matrices 
Ka£inan_ _tor Guin__Matrix 

COL 1 23456 

ROW 

1 9. 11510-01 1,80970+00 -6,86030-01 6.76830-01 -4.57740-01 -7.82920+00 

2 6,32760-03 3.33600-01 5,60230-02 1,66450-02 8,93970-01 3,66190-01 

3 -2,39870-03 5,60230-02 4,57550-01 -2,64430-03 -1,18620-02 6,84870-01 

4 6,76830-01 4.76040+00 -7,56270-01 2,67790+00 2.48630+00 -5,77220+01 

5 -3,84120-04 2,14550-01 -2,84690-03 2. 08640-03 4,71150+00 7,42670-03 

6 -6,57000-03 8,78860-02 1.64370-01 -4,84380-02 7,42670-03 1,87150+00 
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Fi.'udbci£k_G£in Ma£r£x_G^ 


ROW 

COL 1 

2 

3 

4 

S 

6 

1 

4,36130-02 

- 1 .04 17D-01 

-3,35510-01 

3,77920-02 

-8, 18500-02 

-1,38930-01 

2 

0,78490-02 

4 ,04690't'OO 

-9.76770-02 

1 .04270-01 

1 , 28900+00 

-4,60820-02 

3 

-9, 1 1540-04 

6.90760-02 

B. 79310-01 

1 . 35880-04 

3.20020-03 

4.50690-01 


3 Aircraft Poles 

Openj^loo^ £p^es_: X (A ) 

A 



REAL PART 

IMAQ PART 

NAT FREO(HZ) 

ZETA 

FREO(HZ) 

1 

5.3240-01 

0,0 

8. 4730-02 

1 , 000000 

0.0 

2 

3 

1.4250-01 
1 .4250-01 

3 , 7880-01 
-3,7880-01 

6,44 10-02 
6,4410-02 

-0.352163 

-0.352163 

6 .0290-02 
6.0290-02 

4 

6,8200-02 

0,0 

1,0850-02 

1 , 000000 

0,0 

5 

0,0 

0,0 

0,0 

0,0 

0.0 

6 

0.0 

0.0 

0.0 

0,0 

0.0 

^ losed 

-loop_jpoles; 






REAL PART 

IMAG PART 

NAT FREO(HZ) 

ZETA 

FREO(HZ) 

1 

2 

-2,7920+00 

-2,7920+00 

2.8030+00 
- 2 , 8030+00 

6.2960-01 

6.2960-01 

0.705716 

0,705716 

4,4600-01 

4,4600-01 

3 

4 

-1 .8690+00 
- 1 .8690+00 

1 , 8630+00 
- 1 ,8630+00 

4 , 2000-01 
4 , 2000-01 

0,708310 

0,708310 

2.9650- 01 

2.9650- 01 

6 

6 

- 1 , 0640 + 00 

- 1 ,0640+00 

1 ,0360+00 
- 1 .0360 + 00 

2,3630-01 
2 , 3630-01 

0, 7 16447 
0,716447 

1.6490- 01 

1 .6490- 01 
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REAL PART 

IMAG PART 

NAT EHEO(HZ) 

ZEtA 

PREOIHZ ) 

1 

-1.08 10*01 

0.0 

1.7200*00 

1 , OOOOOO 

0.0 

2 

-5.368D+OC 

0,0 

6,5440-01 

1 , OOOOOO 

0,0 

3 

-2.2310*00 

0.0 

3.6500-01 

1 , OOOOOO 

0,0 

4 

-1,1 160*00 

0,0 

1 , 7800-01 

1 , OOOOOO 

0,0 

5 

-6.0100-01 

0,0 

9, 5660-02 

1 . OOOOOO 

0,0 

6 

-4 , 9200-01 

0,0 

7, >1300-02 

1 . OOOOOO 

0,0 


Compensator ^oles 


A A 




REAL PART 

IMAG PART 

NAT FREO(HZ) 

ZETA 

FRE0(H7> 

1 

- 1 ,6790*01 

0.0 

2.6720*00 

1 , OOOOOO 

0,0 

2 

-8,8750+00 

0.0 

1 ,4 130+00 

1 . OOOOOO 

0.0 

3 

-3, 7440*00 

0,0 

5,9590-01 

1 . OOOOOO 

0.0 

4 

- 1 , 1650+00 

0,0 

1 , 8540-01 

1 , OOOOOO 

0,0 

5 

-6,0610-01 

0,0 

9,6460-02 

1 . OOOOOO 

0.0 

6 

-5 , 7070-01 

0.0 

9,0830-02 

1 , OOOOOO 

0,0 


the compensator is stable. 
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